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Effects  of  Terrain  Power  Spectral  Density  Shaping 
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I .  Robert  Ehrl Ich 
Hay  1972 

A  two-dimensional,  f ive-degree-of-f reedom,  digital  computerized, 
wheeled-vehicle-rlde  simulation  is  tested  for  sensitivity  to  two 
parameters:  the  power  spectral  density  slope  of  computer-generated, 
random  terrain  profiles  and  the  spacing  of  the  discrete  profile  points. 
The  vehicle  ride  simulation  Is  exorcised  over  six  terrain  profiles  of 
different  PSD  slopes  but  Identical  RMS  elevations.  The  simulation  Is 
also  exercised  several  times  over  one  basic  profile  described  by  samples 
taken  at  different  measurement  intervals.  Calculated  absorbed  power 
at  the  vehicle  center  of  gravity  is  used  to  compare  ride  roughness  over 
the  different  profiles.  The  vehicle  simulation  Is  shown  to  be  extremely 
sensitive  to  changes  In  PSD  slope.  The  sensitivity  to  changes  In 
measurement  Interval  Is  shown  to  be  dependent  on  vehicle  size  and  mass. 
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INTRODUCTION 


For  some  time,  off-road  mobility  has  bean  of  manlfast  Importance 
to  designers  of  military  vehicles.  The  problem  has  bean  to  design  a 
vehicle  to  traverse  terrain  with  certain  characteristics  at  a 
reasonable  speed  with  negligible  adverse  effects  on  vehicle  payload 
or  driver. 

This  study  concerns  Itself  with  soma  of  the  problem*  encountered 
In.  describing  and  measuring  terrain  characteristics  —  In  particular, 
ground  roughness. 

BACKGROUND 

The  characteristics  of  a  terrain  can  be  numerous.  They  may 

include,  among  others,  slope,  obstacles,  vegetation,  soli  strength, 

and  roughness.  It  has  been  shown,  however,  that  the  single  most 

spend- I Iml ting  aspect  of  off-road  mobility  is  ride  dynamics.*  The 

study  of  ride  dynamics  concerns  itself  with  human  and  cargo  response 

to  vibration.  This  vibration  may  be  caused  by  terrain  roughness 

and/or  discrete  obstacles,  and  filtered  by  the  vehicle  mass,  geometry 

and  suspension  system.  Vibration  caused  by  traversal  of  discrete 

obstacles  (logs,  small  ditches)  Is  transient  in  nature  and  vehicle 

speed  Is  limited  more  by  vehicle  strength  and  operator  Judgement  and 
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experience  than  by  vibrational  characteristics.  On  the  other  hand, 
vibration  due  to  stationary  ground  roughness  is  close  to  a  steady-state 
condition;  dependent  on  vehicle  velocity.  An  operator  will,  if 
properly  motivated,  increase  the  vehicle  speed  until  some  degree  of 


discomfort  Is  felt.  Ho  will  thon  docrosso  speed  slightly.  The  maxl- 
mum  speed  e  vehicle  with  driver  cen  maintain  over  a  certain  terrain, 
then,  Is  determined  by: 

(1)  Terrain  roughness 

(2)  Vehicle  strength  and  suspension 

(3)  Driver  discomfort  (or  payload  delicacy) 

It  has  long  been  recognized  that  seme  quantitative  measure  of  the 
last  of  these,  driver  discomfort,  Is  necessary  to  ride  dynamics  re* 
search.  Many  studies  have  been  conducted  In  an  attempt  to  determine 
some  subjective  measure.  Van  Deusen  shows  the  apparent  futility  of 
this  approach  and  some  of  the  resultant  confusion  In  a  composite 
graph,  reproduced  In  Figure  I.2 

Several  quantitative  measures  have  been  used  with  greeter  success. 
They  Include,  but  are  not  limited  to:  RMS  acceleration  at  the  driver's 
seat  (or  the  area  under  the  RMS  acceleration  time  history)*;  maximum 

p 

acceleration  at  the  driver  seat  ;  the  amplitude  probability  distri¬ 
bution  of  driver  acceleration*;  end  absorbed  power,  a  concept  generated 
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by  Lee  end  Pradko  which  is  Intended  to  quantify  the  energy  dissipated 

by  the  human  In  the  vibration  of  his  limbs  and  flesh  and  counteracted 

partially  by  muscular  control.  In  the  words  of  Lee  and  Pradko  : 

The  Important  characteristics  of  absorbed  power  are 
are  that  it  has  physical  significance  and  therefore 
cen  be  measured  as  well  as  computed  analytically;  and 
that  since  power  is  a  scalar  quantity,  absorbed  power 
can  be  summed  In  complex  multidegree  of  freedom  systems 
to  determine  human  response. 

The  results  of  an  earlier  effort  by  the  same  authors  show  that 
absorbed  power  levels  measured  at  the  driver's  seat  seem  to  agree  more 


1.  Zlegenruecker  and  Mag Id  "Actual  bodily  harm 

fearad" 

2.  Magld  t  Coermann  "One  minute  tolerance" 

3.  Dlekmann  "Intolerable" 

4.  Gorrlll  &  Snyder  "Alarming" 

5.  Parks  and  Snyder  "Alarming"  (3/4"  felt  pad) 

6.  Zeller  "Upper  tolerance  limit" 

7.  Goldman  "Intolerable" 
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Figure  I.  Subjective  Ride  Criterion 
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closely  with  subjective  ride  evaluations  then  do  RMS  accelerations. 

Thus,  the  current  study  will  use  the  absorbed  power  concept.  Even 
among  proponents  of  this  criteria,  there  Is  confusion  over  the 
question  of  what  level  of  absorbed  power  should  be  considered  speed- 
limiting.  Six  watts,  being  a  generally  accepted  figure,  will  be  used 
herein  with  no  apologies  to  those  who  consider  it  too  low. 

This  research  concerns  Itself,  not  with  the  driver  discomfort 
or  vehicle  systems  problem,  but  with  the  first  factor  mentioned  above 
—  terrain  roughness. 

Recently,  (1972)  power  spectral  density  (PSD)  estimates  have  been  used 
to  describe  terrain  roughness  statistically.  The  mathematical 
definitions  of  PSD  will  be  discussed  in  some  detail  later.  When  PSD 
estimates  for  real  terrain  waveforms  are  required,  the  problems  of 
aliasing  and  stationarity  become  important.  Alias!,  g  is  a  condition 
where  false  evidence  of  certain  frequencies  appears  because  the 
discrete  sampling  interval  is  longer  than  the  shortest  wavelength 
present.  Stationarity,  briefly,  is  the  absence  of  long-range  changes 
In  underlying  statistical  properties. 

The  first  proposal  of  the  PSD's  use  to  characterize  stable  ground 
roughness  came  from  Kozin,  Cote,  and  Bogdanoff  in  1963.^  They  found 
that  for  "visually  constant"  ground  roughness,  stationarity  could  be 
safely  assumed  for  terrain  segments  up  to  2000  feet.  They  used 
running  averages  to  correct  for  unavoidable  long-range  trends  such  as 


Others  have  shown  that  the  PSD  curves  of  most  natural  and  man-made 
surfaces  can  be  approximated  with  the  equation: 

PSD(f)  -  Cfn  (I) 

where  f  -  frequency 

which  defines  a  straight  line  on  a  log- log  plot,  with  slope  of  n. 

6 

Van  Deusen,  In  J967»  gave  evidence  to  show  that  for  most  natural  or 
man-made  surfaces,  the  slope  Is  roughly  -2. 

Other  researchers  have  subsequently  used  the  -2  approximation  In 
ride  dynamics  work,  using  computer  simulation  of  both  the  terrain  and 
the  vehicle.  Before  substantive  work  using  a  computer  simulation  Is 
undertaken,  the  assumptions  made  when  Implementing  that  simulation 
must  be  shown  not  to  significantly  affect  the  results  of  the  research. 
In  past  studies,  several  assumptions  have  been  made  regarding  the 
computer-generated  terrain  profiles: 

(1)  That  the  profiles  pass  the  test  for  stationarlty. 

(Murphy  used  a  program  called  STANOR  to  actually 
test  his  prof  I les.) * 

(2)  That,  as  mentioned  above,  the  slope  of  the  PSD 
curve  Is  -2.0  for  any  real  terrain. 

(3)  That  the  input  spacing  can  be  set  to  any 
arbitrary  val'je  (Murphy  used  3.07  inches; 

Kozin  used  up  to  2  feet  to  estimate  the  PSD's 
of  rtai  terrain^). 


The  first  of  these  assumptions  seems  reasonable;  especially  If  the 
input  Is  generated  from  random  numbers  which  are  stationary.  This 
will  be  discussed  In  more  detail  later.  The  second  two  leave  some 
questions  unanswered. 

In  recent  work,  as  yet  unpublished,  Murphy  has  measured  over 
sixty  actual  ground  roughness  profiles.  He  has  found  that  the  slopes 
of  their  PSD  curves  vary  from  -0.6  to  -2.3.  Does  the  assumption  of  a 
-2.0  slope,  then,  cause  significant  errors  in  ride  analysis?  To 
answer  this  question  Is  the  first  objective  of  this  research.  This 
objective  will  be  attacked  by  testing  the  vehicle  ride  simulation  for 
sensitivity  to  changes  in  PSD  slope,  keeping  other  statistical 
parameters  constant. 

The  second  assumption  concerns  a  more  practical  matter.  An 
Investigator  Involved  in  taking  field  data  for  terrain  roughness 
analysis  has  to  trade  off  time  and  money  considerations  against 
accuracy.  Taking  data  points  using  transit  and  rod  Is  time  consuming 
»  thus  expensive.  Taking  data  closer  together  Increases  cost;  while 
choosing  too  great  a  measurement  Interval  decreases  accuracy. 
Theoretically,  some  point  exists  where  further  reduction  of  interval 
wastes  money.  Likewise,  at  some  point  a  further  Increase  causes  a 
serious  reduction  in  accuracy.  The  second  objective  of  this  report 
is  to  obtain  some  quantitative  grasp  of  the  accuracy  of  vehicle  model 
performance  versus  distance  between  data  points.  This  will  be 
accomplished  by  creating  a  single  terrain  profile  with  a  small 
Interval  and  exercising  a  computer  vehicle  simulation  over  this 
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DISCUSSION  AND  COMPUTER  CALCULATIONS 


NOIPSD 

Th«  first  requirement  of  this  study  was  the  generation  of  terrain 
profiles,  created  with  random  amplitudes.  In  addition,  some  control 
of  statistical  properties  was  necessary.  A  computer  program  (NOISE  1 ) 
written  for  this  purpose  was  obtained  from  Mr.  N.  S.  Murphy  Jr.  of  the 
Waterways  Experiment  Station.  Another  program,  PSD,  was  designed  to 
accomplish  the  PSD  estimation.  To  reduce  the  execution  time,  NOISEI 
and  PSD  were  combined  to  form  the  program  NOIPSD.  Unfortunately,  both 
were  written  for  a  GE400  machine,  while  the  computer  available  was  a 
PDP-10.  This  necessitated  several  adjustments  In  the  program.  In 
addition  to  translation  Into  a  different  Fortran  language.  A 
description  of  the  program  and  the  changes  made  In  Its  basic  structure 
follows.  A  program  Using  Is  In  Appendix  A. 

First,  a  random  number  generator  was  required.  Since  these 
generators  (sub* prog rams)  are  generally  machine-specific,  some 
alteration  of  the  program  was  necessary.  An  Internal  PDP-10  program, 
called  RAN(Z),  was  used  to  obtain  twelve  uniform  random  numbers 
between  zero  and  one  for  each  profile  point.  It  was  desired  that  each 
terrain  profile  be  different.  RAN(Z),  however,  has  a  fixed  starting 
number  for  each  program  execution.  This  Is  analogous  to  having  a 
fixed  table  of  random  numbers  from  which  to  choose.  Obviously,  then, 
each  execution  would  result  In  the  same  profile  shape.  It  was  decided 
to  use  a  scheme  In  which  a  "skipping  number"  was  used  to  cause  the 
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program  to  produce  different  profiles  while  melntalnlng  their  random 
nature.  IWelve  normal  random  numbers  between  zero  and  one  ere  used 
to  obtain  random  numbers  with  zero  mean  by  the  uniform  deviates 
method: 

Vj  -  on  (1%,  -  6.0)  (J  -  1,2,3 . 1200) 

where  Vj  ■  a  random  normal  number  with  zero  mean 

Zj  ■  a  uniform  random  number  from  RAN  (Z) 

a  »  desired  standard  deviation  of  V 
n 

Before  each  twelve  numbers  used  from  the  "table"  of  random  numbers, 

N6  (the  "skipping  number")  points  were  skipped.  N6  should  be  kept 
relatively  low  to  reduce  computer  time.  Prime  numbers  were  assigned 
to  N6  each  time  N0IPSD  was  run  to  Insure  profile  difference. 

After  all  the  Vj's  are  determined,  their  mean  Is  computed.  It 
should  be  zero  and  usually  is  very  small.  To  make  sure,  each  Vj  is 
then  shifted  by  the  computed  mean. 

This  procedure  gives  white  Gaussian  noise,  which  has  essentially 
a  level  power  spectral  density  curve.  The  above  spectrum  Is  shaped 
using  a  digital  simulation  of  an  analog  low-pass  filter: 


where  Y|  •  the  resulting  profile 

a  ■  the  special  cutoff  frequency 
t  >  the  interval  between  points  (constant) 
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Th«  product  err  Is  usod  to  adjust  ths  frequency  content  of  the 
prof I  Is.  For  most  of  this  study  ,  r  was  kept  at  a  constant  4.0 
Inches.  The  value  chosen  for  or  ,  then,  determined  the  ultimate 
power  spectral  density  slope.  A  table  of  or  values  was  determined 
by  trial  and  error  and  Is  given  below: 

Table  I 

(Valid  for  t  -  4.0) 


Oes 1  red 

Use  NG 

Use  a 

Slope 

of 

of 

-0.60 

2 

.2575 

•1.20 

19 

.0888 

-1.85 

13 

.0129 

-2.00 

7 

.0052 

-2.15 

II 

.00138 

-2.30 

3 

.00016 

Next,  the  desired  RMS  level  (DRMS)*  is  achieved  by  computing  th* 
actual  RMS  (ARMS)  and  adjusting  each  point  by  the  factor  DMS/ARMS. 
This  results  In  the  desired  results,  and  completes  the  profile  gen¬ 
eration  portion  of  the  program. 

Before  continuing  on  In  the  program  description,  the  mathematical 
definitions  of  the  autocovariance  function  and  power  spectral  density 
estimation  need  to  be  presented. 


* 


Names  In  parentheses  indicate  program  variables.  See  Appendix  A. 
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Assume  X  ( t )  Is  a  continuous.  Inf  In  Its  waveform.  The  covarlancss 
(Cjj)  art  defined  by • ^ 

C,j  -  cov  |x(t,),X(tj)| 

-  ave  |lX(tJ)-Sr(t|)]*lX(tj)-5T(tj)j| 


where  X  Indicates  the  mean. 


Since  the  terrain  waveforms  are  created  with  zero  mean,  3f(tj)  - 
X(tj)  -  X(tR)  ■  0.0.  It  follows  that: 

cu  -«•{*<«, 

Furthermore,  assuming  statlonarlty,  the  covariances  will  depend  on 


time  separation  only. 


C,j  -  C(t,-tj)  -  C(t) 


whore  T  Is  generally  called  lag.  From  this,  and  assuming  ergodic 
properties,  the  covariance  at  lag  t  can  be  written: 


C(t)  ■  ave  jx(t)*X(t+T)| 


or.  In  functional  notation: 


■  "■{  r 

T—  J., 


X(t)*X(t+T)dt 


which  Is  generally  called  the  autocovariance  function.  Equation  (2) 
may,  In  the  Ideal  case,  be  reduced  to:^ 


ft- 106 


Sine*  C(t)  If  an  aven  function  (symmatric  about  tha  zaro  axis),  the 
aquation  may  ba  simply  written: 


In  tha  particular  casa  cons  I  da  rad  herein,  tha  data  to  be 
analyzed  Is  discrete  and  equally  spaced.  One  very  Important  problem 
to  be  '*nnsidered  *n  these  cases  is  aliasing.  If  the  data  Is  taken 
from  a  continuous  waveform  at  equally  spaced  intervals,  as  would  ba 
the  case  In  recording  a  real  terrain,  Some  frequency  estimates  are 
distorted.  The  basic  problem  Is  illustrated  in  Figure  2  which  shows 


ft- 106 


\k 


how.  In  equl-spaced  sampling  of  a  stna  wava,  an  Imaginary  longar 
wavalangth  appaars  to  ba  prasant.  In  tha  prasant  cata,  tha  dlscrata 
tarraln  data  ara  not  naasurad  from  a  continuous  sourca,  but  craatad 
in  Its  dlscrata  form.  Tha  highest  frequency  possible  Is  l/(2At)  which 
Is  usad  as  tha  high  cut-off  frequency.  Aliasing,  In  this  special 
case,  Is  of  no  concern. 

For  discrete  data,  Blackman  and  Tukay  suggest  a  thraa-stap  method 
for  power  spectral  density  estimation:^ 

(1)  If  X(J),X(2),X(3),...,X(N)  Is  the  discrete  series  with 
Interval  A1  between  values,  tha  autocovariance  Is 
computed  with  lag  t  ■  h*Ai  as  follows: 

C(T)  ■  n=?7F  q Jfo  (x(q)*x(q-^r>)  (5) 

where  N  -  number  of  points 
r  —  0, i,2,««.,m 
m  <  N/h  ■  maximum  lag 
h  >  0  (integer) 

(2)  The  raw  spectral  estimate  Is  computed: 

m-l  / 

Mc 

where  f(r)  ■  r/(2mAT) 

r  -  0, t,2,...,m 
C,  -  C(Tf) 


*cos(qrTT/m) 


♦  Cm*cos(rn) 


>1 


J 


(6) 


V(r)  -  At  C  +2 
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(3)  Du«  to  tho  flnlta  n unbar  of  data  to  ba  analyzad  and  tha 

disc rata  nunbar  of  lags,  tha  spactral  astlmata  nust  ba  ' 

snoot  had.  Hanning  Is  usad  In  this  analysts: 


u<f>  ■  Al.*r  ♦  jll  A.J  [vr*J  ♦  Vj]  »> 

nhara,  for  banning:  (l«3) 

^30  " 

Aj,  -  0.23 

Ajj  •  0.00  for  J  -  2,3,4,... 


to  turning  to  tha  computar  program,  N0IPSD,  tha  first  stap  In  tha 
PSD  astlmatlon  procass  Is  to  conputa  tha  maximum  nunbar  of  lags  (LLAG). 
From  Blackman  and  Tukay  tha  aquation  ls:^ 

K  »(f  •  *$*) 

LLAG 

wham  K  ■  tha  aqulvalant  n unbar  of  dagraas  of  fraadom  of  a  chl-squara 
distribution 

N  ■  nunbar  of  points  In  prof I  la 
LLAG  ■  nunbar  of  lags  (max  lags) 
v  ■  prof I  la  spacing 

Solving  tha  abova  aquation  for  tha  numbtr  of  lags: 

/ 


....  6n 
IWS  *  TjJiiJiy 


R- 161*6 


16 


Obviously,  a  balance  axlats  among  tha  number  of  points,  N,  the 
number  of  lags,  and  K.  The  choice  of  K  determines  the  degree  of 
accuracy  of  tha  estimate,  as  seen  tn  Figure  3.  It  shows  the  distri¬ 
bution  of  FSD  estimates,  for  Instance,  as  fixed  multiples  of  their 
average  values.  As  an  example,  consider  a  PSD  estimate  which  has  an 
average  value  of  10  InVcycle.  For  K  -  10,  Individual  estimates  will, 

In  the  long  run,  fall  below  .49  times  Its  average  value  (4,9  I n^ /cycle). 
In  brief,  80  percent  of  all  values  would  fall  within  the  Intervrl  4.9 
to  16.0.  As  can  be  seen,  higher  values  of  K  give  higher  accuracy.  But 
for  a  fixed  N,  the  numtx  r  of  lags  decreases  with  K,  reducing  the 
amount  of  Information  In  the  PSD  curve.  Thus  with  a  very  high  K,  one 
point  on  the  PSD  curve  would  result.  It  would  be  vary  accurate,  but 
not  of  much  use.  For  N  -  1200  and  t  -  4.0,  a  value  of  K  -  20  gives 
30  lags.  This  value  of  K  gives  reasonable  accuracy  of  estimates  while 
giving  a  reasonable  number  of  lags.  Since  the  program  must  allow  for 
changes  In  N  and  t  ,  the  equation  was  written: 

m  -  LLAG  -  6N/(62*  t) 
which  Is  used  In  the  program. 

The  autocovariance  computation  is  done  by  using  h»1  In  equation 
(5).  A  check  Is  made  on  the  RMS  level  at  this  point  by  taking  the 
square  root  of  the  first  (r«0)  autocovariance,  which  should  equal  the 
desired  RMS. 

Next,  the  raw  spectral  estimates  (PX)  are  computed  using 
aquation  (4),  and  smoothed  by  hamming  to  give  the  smoothed  power 
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spectral  density  estimates  (SPX): 

SPX(I)  -  .54PX(1)  ♦  ,46PX(2) 

SPX(J)  -  .23PX(J-I)  ♦  .S4PX(J)  +  .23PX(J+1) 

(J“2,3»^»»»»  ,rn-0 

SPX(m)  -  .54PX(m)  +  ,46PX(m-l) 

Finally,  after  discarding  PSD  points  outside  the  low  and  high 
cut-off  frequencies,  the  smoothed  PSD  estimates  are  fitted  to  a 
straight  line  rn  a  log-log  plot,  using  a  least  square  fit  routine. 

The  result  Is  an  equation: 

PSD  -  Cfn 

where  C  -  the  Intercept  of  the  line  at  f-1.0 
f  ■  frequency 
n  -  slope  of  the  line 

This  completes  the  discussion  of  NOIPSD. 

VEH 

* 

The  next  major  requirement  of  this  work  was  a  vehicle  simulation. 
Again  this  was.  In  Its  original  form,  obtained  from  Hr.  N.  S.  Hurphy 
Jr.,  WES.  The  program  was  also  translated  to  POP-IO  FORTRAN  IV  for 
use  herein.  Several  subprograms  were  changed  or  added  for  this  study, 
Including  DATA,  GAMSUB,  WHEELS,  PRINT,  and  FILIN.  Basically,  this 
computar  program  Is  a  2-dlmens Iona),  5-degree-of-freedom  simulation. 
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Th«  following  assumptions  war*  made  In  deriving  the  equations  of 
motion: 

(1)  The  vehicle  sprung  mass  Is  rigid. 

(2)  The  only  forces  acting  on  the  sprung  mass  are 
suspension  forces. 

(3)  There  Is  no  surge  acceleration  or  motion  In  sway, 
yaw,  or  roll  degrees  of  freedom. 

(4)  Pitch  Is  small  enough  to  allow  a  small-angle  approximation. 

Assumption  3  was  made  reluctantly,  but  was  necessitated  by  time 
and  cost  considerations.  The  others  are  reasonable  considering  the 
purpose  of  the  simulation. 

Figure  4  shows  a  schematic  of  the  general  wheeled  vehicle  model. 
Three  axles  are  shown  and  that  is  the  maximum  number  which  can  be 
simulated  with  this  model.  A  two-axle  vehicle  Is  modeled  by  setting 
Jtj,  and  mj  equal  to  zero.  Thus,  the  one  model  can  simulate  almost 
any  conventional  vehicle.  In  this  study,  ont'  2-axle  and  one  3-axle 
vehicle  are  used.  They  are  the  HIS  I  Jeep  and  the  M33  2$T  Truck, 
respectively. 

The  equations  of  motion  were  written  using  Newton's  Law  on  each 
of  the  free-body  masses  shewn  In  Figure  5. 

Body: 

Sum  of  vertical  forces: 

V  "  Fl  +  F2  '  V 


(8) 


tin  model 


Schematic  of  General 
Wheeled  Vehicle  Model 
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Sum  of  moments  about  pitch  axis: 

ie  -  F,i,  co»  e  -  f 2'i2  cos  e  (9) 

Unsprung  masses: 

Sum  of  vortical  forces: 

mjk'j  -  -Fj  -m,g  +  Pj  (1-1,2, 3)  (10) 


The  angle  0  Is  defined  In  the  rear  suspension  geometry  as  In 
Figure  5.  The  forces  F^F^  •  ar><*  Pj  are  c'^f Ined: 

F,  -  k,*,-*,^ 

A|  -  ij-z-Jfcj  sin  6 
A|  -  Zj-Z-ij  6  cos  0 


*2 

*2 


Zg-s+^z  »ln  0  -f'j  sin  0 


*2’*+J*2®  cos  0  -^3  0  cos  0 


f}  ■  kjiycjij 

■  Z3"Z+J*2  s,n  ®  +  ^4  ,,n  P 
-  ij-z+^z0  cos  0  ♦  0  cos  0 


Pj  -  forces  from  terrain  profile  as  transmitted  by 
the  tire  model. 
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The  program  named  VEH  Is  designed  to  solve  equations  8,  9,  and 
10.  Basically,  It  performs  Its  function  using  a  fourth  order 
Runge-Kutta-Gt II  algorithm.  A  program  listing  Is  In  Appendix  B. 

Before  the  calculations  are  begun,  program  options  and  variable 
parameters  are  fed  Into  the  program,  normally  by  teletype  on  a  remote 
hookup.  The  program  options  determine  which  parts  of  the  pro^  -m 
will  be  executed.  The  variable  parameters  include  the  Input  profile 
name,  tire  deflections,  vehicle  velocity  In  miles  per  hour  and  the 
teletype  printout  time  Interval.  These  parameters  will  be  explained 
more  fully  In  the  text  that  follows. 

First,  VEH  calls  subroutine  FILIN.  On  this  initial  call,  FILIN 
requests  an  Input  of  des  I  red  A l  (OELTAL)  ,  which  may  be  any  Integer 
multiple  (*1)  of  the  input  profile  spacing.  FILIN  then  reads  the 
actual  input  profile  spacing  (SPACING),  the  profile  Identification 
line  (FIO)  and  the  first  ten  profile  points  (FYIN)  from  the  input 
file.  A  parameter,  HM,  is  then  calculated  by  the  following  formula: 

^  DELTAL 
"  SPACING 

which  defines  the  number  of  profile  points  (FYIN)  to  be  skipped  between 
each  two  points  returned  to  the  main  program.  Thus,  If  the  spacing 
In  the  Input  profile  Is  1.0  Inches  and  the  desired  spacing  Is  4.0 
Inches,  MM  ■  4.  Every  fourth  point  will  be  returned  to  the  main 
program.  This  feature  was  added  to  allow  investigation  In  this  study 

The  words  In  parentheses  are  the  parameter  names  in  the  program 
listing.  See  Appendix  B. 
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of  the  effect  of  Input  profile  spacing  on  simulation  performance. 

After  the  vehicle  name  and  the  program  options  are  read,  the 
program  calls  subroutine  DATA.  DATA  sets  the  vehicle  parameters  to 
the  proper  values.  It  must  contain  every  vehicle  parameter  necessary 
to  simulate  the  desired  vehicle.  Some  of  the  data  statements  contain 
values  of  Interest  only  for  tracked  vehicles.  The  tracked  vehicle 
subroutines  are  not  Included,  but  the  data  was  left  In  case  later 
Investigators  wish  to  use  the  program.  The  tracked  vehicle  subroutines 
are  available  If  desired.  For  wheeled  vehicle  simulation,  a  list  of 
necessary  vehicle  parameters  for  inclusion  In  subroutine  DATA  Is  given 
In  Appendix  C. 

After  most  of  the  parameters  are  set,  but  before  returning  to  the 
main  program,  DATA  calls  subroutine  GAMSUB,  which  computes  the 
parameters  used  in  the  radlal-sprtng  tire  model.  The  tire  model  Is 
constructed  under  the  assumption  that  a  tire  can  be  simulated  by  a 
series  of  radial  springs,  as  in  Figure  6.  In  the  original  program, 
the  calculations  done  in  GAMSUB  were  made  by  hand  and  Inserted  through 
subroutine  DATA.  In  this  study  it  was  desired  to  change  the  Input 
profile  spacing  from  one  execution  to  the  next.  Thus,  GAMSUB  takes 
Into  account  the  input  profile  spacing  (DELTAL)  and  "creates"  radial 
springs  such  that  the  distance  between  the  projections  of  their  outer 
ends  on  a  horizontal  surface  will  be  exactly  equal  to  the  profile 
spacing.  Thus,  Dj,  are  mu*t,P,es  DELTAL.  The  maximum 

angle  allowed  from  the  vertical  is  roughly  S3  degrees.  One  spring  Is 
always  placed  vertical,  the  others  symmetrical  to  front  and  rear. 
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Basically,  Y  (GAMMA)  Is  tha  vertical  component  of  tha  radial 
spring  force-deflection  function.  Each  radial  spring  Is  assumed  to 
have  the  same  spring  rata,  which  Is  determined  separately  for  each 
tire  from  the  tire  force-def lection  curve.  The  curves  for  MI5I  tires, 
for  Instance,  ere  shown  In  Figure  7.  The  curves  are  very  close  to 
being  straight  tines,  so  a  linear  approximation  Is  used.  The  load  on 
each  tire  Is  computed  (note  that  for  3-axle  vehicles  the  rear  tires 
are  assumed  to  be  duals)  and  typed  out,  In  order  from  front  to  rear, 
on  the  teletype.  After  each  load,  a  tire  deflection  Is  entered  by  the 
operator.  Care  must  be  taken  to  enter  tha  deflection  at  the  propef 
Inflation  prassure.  The  M15I  Jeep,  for  Instance,  should  have  18  psl 
In  the  front  tires  and  22  psl  In  the  '  tar  for  cross-country  operation.® 
GAMSUB  then  computes  the  radial  spring  force  for  each  tire.  SPKF, 

SPKRI,  AND  SPKR2  are  the  program  variables  for  these  spring  forces 
front  to  rear.  They  are  computed  In  the  following  manner: 

Suppose  the  dead  load  on  the  tire  causes  a  tire  deflection  (Y) 
as  In  Figure  8.  The  radial  spring  oonstant  (SPK)  Is  computed  from  the 
dead  load  (WEIGHT)  and  spring  daflectlons  (DELTA)  as  follows: 

NSEGS 

WEIGHT  -  £  SPK*cos  <^*DELTA,  (II) 

where  SPK*cos  6j  ■  the  vertical  component  of  SPK 

NSEQS  -  2*KK+1  ,  the  number  of  radial  springs  In  the  tire  model. 

SPK  -  MEGS  W£lf'T -  <”> 

£  cos  tp.^DELTA, 

1-1  1  1 


.-«■  - -  .  i.r,  .  .  .  ,  -  - 


sine*  SPK  is  assumed  constant*  In  addition*  since  cos  6,  can  be 
written  (sea  Figure  6): 


cos  6 1 


equation  (9)  becomes: 


SPK  -  WEIGHT 

NSEGS  D{ 

E  -r*  DELTA 
l-l  * 


03) 


where  DELTA |  -  Y  -  THRESH,  -  z 

(If  DELTA,  Is  <  0.0,  DELTA,  -  0.) 
THRESH,  »  R  -  0'  (see  Figure  9) 

R  >  the  undeflectad  tire  radius 


Note  that  In  this  case,  Y  causes  deflections  only  In  three  of  the 
radial  springs.  Thu*,  DELTA,  -  CELTA2  -  DELTAfi  -  DELTA?  -  0.0  and 
equation  (II)  becomes: 


SPK  - 


WEIGHT 


°j  D4  ^ 

■f  *  DELT^j  +  -J  *  DELTA,.  4  *  DELTA, 


R 


The  radial  spr’ng  rates  for  different  tires  on  the  same  vehicle 
may  vary,  either  because  of  different  inflation  pressures  or  different 
load  distribution. 

As  the  final  step,  GAMSUB  computes  the  values  for  GAMMA  by  the 


equation: 


GAMHAjj  -  SPKj  *  CO*  9, 

1-1. 2, 3 . MSECS 

J-t.2,3  (axle  number) 


For  computer  notation,  the  GAMMA'*  ere  numbered  from  front  to  rear  of 
the  vehicle  as  shown  In  Figure  10* 


Figure  10:  Numbering  Scheme  for  Tire  Segments 


Upon  return  to  the  main  program,  the  number  of  steps  to  be  used 
In  the  Runge-Kutta-GIII  (RKG)  algorithm  Is  computed,  based  of  vehicle 
velocity.  The  step  size  Is  desired  to  be  roughly  .001,  which  Is 
adjusted  slightly  to  Insure  an  exact  number  of  steps  between  Input 
terrain  profile  points.  This  completes  the  preliminary  calculations. 

The  program  enters  an  integration  loop  at  this  point,  starting 
effectively  with  the  calling  for  the  second  time  of  subroutine  FILIN. 
Each  time  FILIN  is  called  from  the  Integration  loop,  it  simply  returns 
one  profile  point  (YIN)  to  the  main  program. 

If  a  detailed  output  file  Is  desired,  the  program  calls  FILWRT, 
which  writes  time  (T),  YIN,  and  absorbed  power  (ABSPWR)  for  each  step. 
It  also  writes  displacement,  velocity,  acceleration,  and  RMS 
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The  fore*  exerted  through  the  suspension  system  to  the  body  must  now 
be  computed  for  each  axle. 

WHEE15  allows  considerable  flexibility  In  modeling  spring- 
deflect  ton-vs. -force  curves  and  In  modeling  damp I ng- force- vs. - 
deflect  I  on- velocity  curves.  Both  curves  are  specified  by  data  state¬ 
ments  In  subroutine  DATA. 

The  spring  force-deflection  function  of  a  typical  suspension 
spring  can  be  approximated  by  five  linear  segments,  as  shown  In 
Figure  II. 

The  deflection  axis  can  be  divided  Into  five  regions: 


Region  1: 

—CD 

to 

X1 

Region  2: 

to 

X2 

Region  3: 

X2 

to 

X3 

Region  4: 

x3 

to 

x4 

Region  5: 

x4 

to 

+CD 

The  X|  ere  the  region  limits  (SLIMIT|). 

In  eech  region,  the  force-deflection  function  Is  approximated  by 
a  linear  equation  of  the  form: 

FORCK,  -  m  *  SPDEF,  +  C 

where  FORCK  ■  rhsultant  spring  force 

m  -  slope  of  the  line  (SSLOPE) 

SPDEF  ■  deflection  of  the  spring 

C  -  Intercept  of  the  line  at  SPDEF  ■  0.0  {SINT) 

|  -  axle  number 


acceleration  of  tha  cantor  of  gravity,  tha  pitch,  and  each  axle.  If 
the  talatypa  printout  Interval  la  exceeded,  subroutine  PRINT  la 
called,  which  causes  the  same  Information  to  be  typed  out.  If  desired 
PRINT  Mill  cause  only  time  and  absorbed  power  to  be  printed  out. 

PRINT  also  asks  the  operator  If  he  desires  to  stop  execution.  If  the 
answer  Is  yes,  the  main  program  transfers  out  of  the  Integration  loop. 
If  tha  answer  Is  not  yes,  the  main  program  then  calls  subroutine  SHIFT 

SHIFT,  as  the  name  Implies,  causes  each  profile  point  to  be 
shifted  by  DELTAL  inches  to  the  rear  of  the  vehicle.  The  main  program 
then  sets  the  first  profile  point  to  equal  YIN  and  continues. 

Next,  the  differential  equations  are  solved.  For  each  M(G  step 
between  profile  points,  the  suspension  forces  are  calculated  by  a 
subroutine  celled  WHEELS. 

The  first  thing  done  by  WHEELS  Is  the  computation  of  the  forces 
on  each  axle  (FONCW).  It  uses  the  same  method  as  before  to  determine 
radial  tire-model  spring  deflections,  OELTA: 

DELTA,  -  Yk  -  THRESH,  -  Zj 

where  J  "axle  number 

no  negative  values  of  DELTA  are  allowed 
Y^  ■  elevation  of  the  profile  point  under  the  k-th  spring 

The  axle  forces  (FORCW)  are: 

NSEGS 

FORCW.  -  £  DELTA  t*GAMMA,  (J-l,2,3) 

J  |-|  11 
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The  simulation  of  r>e  suspantlon  damping  function  Is  assantlalty  tha 
sama.  Data  I  Is  ara  glvan  In  Appendix  C. 

Thus,  aftar  computing  tha  suspanslon  spring  daflactlon  (SPDEF) 
and  daflactlon  velocity  (OSPDf),  subroutine  WHEELS  compares  their 
values  with  tha  region  limits,  decides  which  aquations  to  use,  and 
computes  tha  suspension  forces  (FORCK  and  DAMP). 

Execution  then  returns  to  the  main  program,  which  Immediately 
calls  subroutine  RUNGE,  tha  RKG  Integration  schema.  Next,  If  the 
absorbed  power  option  Is  requested,  subroutines  POWER  and  RUNGE  are 
used  to  compute  It. 

The  concept  of  absorbed  power  in  mobility  research  appeared  as 
an  attempt  to  resolve  the  confusion  In  driver  vibration  limits  ex- 

3  L 

hlblted  by  Figure  I.  It  was  theorized  by  Lee,  Pradko  and  others  ' 
that  the  driver  will  adjust  the  vehicle  conditions  (e.g.,  speed)  so 
that  he  can  completely  absorb,  by  flexing  and  unflexlng  his  muscles, 
all  the  power  of  the  vibrations  he  receives  from  the  vehicle  In  order 
to  keep  his  eyes  or  hands  steady  to  see  clearly  and  operate  the 
vehicle  controls.  A  series  of  experiments  indicated  a  great  deal  of 
merit  In  the  concept  and  determined  that  6  watts  was  a  general  level 
of  power  that  the  human  can  or  is  willing  to  absorb  during  driving. 

This  has  been  accepted  by  many  investigators  In  ride  dynamics,  not 
so  much  as  an  ultimate  truth  but  as  the  best  available  now. 

In  the  course  of  these  studies  several  ways  of  calculating 
absorbed  power  were  developed  (the  first  three  reported  in  Lee  end 
Pradko^): 


1 
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A*  For  Infinite  averaging  time: 

avarago  absorbed  power  ■Urn  *  f  F(t)v(t)dt 

T—  TJ0 

where:  F(t)  -  Input  forca  to  tha  drlvar  support  (e.g.,  taat) 

V(t)  ■  Input  valoclty  of  tha  drlvar  support 

B.  For  a  finite  averaging  tlma: 

—  ^  (f*)  ♦  p(t)  -  Kf(t)  v(t) 

S)  dr  dt  n 
n 

P(t)  ■  finite  avarago  absorbad  power 
A  ■  damping  factor 

«D  •  lowest  frequency  to  be  averaged  (rad/sec) 
n 

F(t)  -  Input  forca  to  tha  drlvar  support 
V(t)  ■  Input  velocity  to  tha  drlvar  support 
K  *  conversion  constant 

C.  In  tha  frequency  domain,  It  can  be  computed: 

N 

P-  E  K(f.)(RMSA(f?)) 

1-0 

where  P  ■  finite  average  absorbed  power 

K(fj)  -  parameter  dependent  on  frequency  I 
RM$A(f|)  ■  RMS  acceleration  at  frequency  I 

Tables  of  K(fj)  are  given  In  Reference  3* 

0.  Recently,  work  has  bean  dona  at  WES  to  create  an  algorithm  to 
compute  absorbed  power  In  tha  time  domain.  Tha  result  was  tha  following 


Mt  of  equations,  derived  from  the  analog  circuit*  shown  In  Flgura  12. 
a(t)  Is  tha  acceleration  of  the  driver  support,  the  other  variables 
say  be  considered  as  Intermediate  quantities  needed  In  tha  calculation 
of  absorbed  power. 

&n  -  -0.1755  a(t)  -  2.742  ^ 

*10  "  “,#755  ‘  38808  x||  “  h6t67  xio 

U2  -  -0.1755  a(t)  -  1.042  x10 

-  -10  Ug  -  6.249  Uj 

kB  -  -10  Ug  -  78.59  ^  -  55.28  Xg 

U1  ■  Ug  -  3.246  Xg 

Kj  -  -100  U]  -  47.78  uo 

ftg  -  -10  Uf  +  71.6  Kj  -  53.49  Xg 
u0  -  -  U,  +  1.318  Xg 
*5  -  -100  u0  -  59  x5 
Xg  -  -0.01294  a(t) 

X|  -  0.00073  Xg  x5 

100  x. 

PWR  -  — 

In  a  study  yet  to  be  published,  Hurphy  has  found  this  latter  algorithm 
to  be  accurate  within  the  frequency  range  of  Interest  (.50  to  20  cps). 
His  version  of  subroutine  POWER  was  used  In  the  program. 


After  the  Integration  schema  Is  completed  for  each  step,  the 
output  scaling  Is  accomplished,  RMS  accelerations  are  calculated, 
and  a  new  profile  point  Is  obtained  from  FILIN  to  begin  the  Inte¬ 
gration  loop  anew. 


multiplier 


Absorbed  Power 
I  volt  »  /  well 


Analog  Ci rcui t  for 
Absorbed  Power  Computation 
given  by  Murphy* 
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TEST  PROCEDURE  AMO  RESULTS 

L  The  first  objective.  as  stated  previously,  was  to  tast  tha 
sansltlvlty  of  tha  vahlela  simulation  to  changes  In  tha  PSD  slopa  of 
tha  Input  tarraln  profits. 

Six  tarraln  profllas  wars  created,  using  MO  I  PSD,  with  PSD  slopes 
ranging  from  -0.6  to  -2.3.  Each  profits  had  1200  points  spacad  at 
4-Inch  Intervals  and  an  RHS  elevation  of  4.0  Inches.  Tha  Inputs  to 
N0IPSD  for  creation  of  these  profiles  are  repeated  In  Table  2. 


Table  2 

(valid  for  t  ■  4.0) 

Desired 

Use  NG 

Use  or 

Slope 

of 

of 

-0.60 

2 

.2575 

-1.20 

19 

.0888 

-1.85 

13 

.0129 

-2.00 

7 

.0052 

-2.15 

II 

.00138 

-2.30 

3 

.00016 

For  PSD  estimation,  a  low  cut-off  frequency  (FLOW)  of  .0002  was  used 
In  each  case.  This  caused  the  zero-frequency  estimation  to  be 
deleted  while  all  others  were  used  In  the  curve-fitting  routine.  The 
zero-frequency  (or  infinite  wavelength)  component  of  the  terrain 
profile  would  have  no  effect  on  the  vehicle  simulation.  The  resultant 
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PSD  data  and  aquations  ara  shown  In  Figures  13  through  18.  Flgura  19 
shout  all  six  aquations  superimposed  for  comparison. 

Flguras  20  through  22  show  tha  Initial  40  faat  of  tha  tarraln 
prof  I  las.  Tha  vortical  sea  las  of  the  profiles  ara  dlstortad  by  a 
factor  of  sight,  so  they  appear  rougher  than  they  should.  The  high 
frequency  component  can  be  seen  to  decrease  as  the  PSD  slope  becomes 
larger,  which  conforms  to  the  expected  trend. 

Originally  each  vehicle  was  to  be  exercised  over  each  profile, 
and  by  trial  and  error,  a  "critical  speed"  determined.  This  critical 
speed  was  defined  as  that  spead  at  which  the  vehicle  exhibited  six 
watts  absorbed  power  at  the  driver's  seat.  The  M151  Jeep  simulation 
was  exercised  first.  Each  simulation  was  run  until  the  absorbed 
power  seemed  to  stabllllse  (until  several  consecutive  values  were 
very  close).  This  normally  took  from  4-6  seconds  of  real  time.  The 
following  results  were  thus  obtained: 

Table  3 


PSD 

M151 

Slope 

Critical  Speed 
(mph) 

-0.60 

2 

-1.20 

2 

-1.85 

4 

-2.00 

5 

-2.15 

10 

-2.30 

17 

This  data  Is  plotted  In  Figure  21.  At  the  end  of  one  of  the  last  runs, 


Section  of  Profiles  1  and  2 
horizontal  scale  in  feet) 


Figure  21.  Section  of  Profiles  3  and  4 
(N01E:  horizontal  scale  in  feet) 


Section  of  Profiles  5  and  6 
horizontal  scale  in  feet) 


however,  the  absorbed  power  Jumped  up  sharply.  To  determine  whet  may 
have  happened,  the  M15 1  was  exercised  over  one  entire  terrain  (slope 
of  -2.3)  at  a  velocity  of  eleven  miles  per  hour.  The  absorbed  power 
vms  plotted  using  a  special  computer  program,  PWRPLT  (see  Appendix  D). 
For  this  execution  the  absorbed  power  leveled  off  at  a  value  of  5.3 
watts  after  about  18  seconds  of  ival  time,  as  compered  to  where  It 
began  to  level  off  at  3.4  watts  around  6  seconds.  Thus,  the  critical 
speeds  determined  for  the  M151  Indicate  the  sensitivity  of  the 
simulation  but  can  be  assumed  to  be  about  50  percent  too  high. 

Since  It  took  four  or  five  runs  to  find  each  critical  speed,  It 
was  determined  that  to  continue  the  "critical  speed"  approach  would 
Increase  the  cost  of  the  program  by  a  factor  of  from  3  to  5.  Since 
the  necessary  funds  were  not  available,  a  slightly  different  approach 
was  taken. 

First,  each  vehicle  was  exercised  over  each  terrain  profile  at 
5  mph  and  at  18  mph.  These  two  speeds  are  near  the  relative  minimum 
and  maximum  of  off-road  mobility  expectations  and  should  exercise 
each  vehicle  In  most  of  Its  tolerance  range. 

Second,  Instead  of  measuring  discrete  absorbed  power  values,  a 
new  subprogram,  AVERAGE  was  written  to  compute  a  running  average  and 
PWRPLT  was  usad  to  plot  the  average  absorbed  power. 

Each  vehicle-profile  combination  underwent  a  15-20  second 
simulation.  The  final  average  absorbed  power  was  used  to  determine 
simulation  sensitivity  to  the  PSD  slope  of  the  Input  terrain  profile. 


The  results  are  given  In  Figures  24  end  25. 

Unfortunately,  at  both  5  mph  end  18  mph  on  the  terrains  with 
PSO  slopes  of  -0.6  and  >1.2,  both  vehicles  exhibited  large  pitch 
verlatlons.  Since  a  small  pitch  angle  approximation  was  used  to 
derive  the  vehicle  simulation,  ell  data  on  those  profiles  was 
considered  Invalid  and  Is  not  Included  In  the  analysis. 

There  Is  no  scientific  reason  to  draw  a  straight  line  between 
the  two  points  at  each  slope.  In  fact,  according  to  Hurphy's  data,* 
the  lines  should  probably  be  concave  downward  (second  derivative  of 
the  power  with  respect  to  slope  negative).  The  lines  are  drawn  simply 
to  enhance  visual  comparison.  Even  a  hasty  Inspection  makes  one  thing 
obvious  --  the  vehicle  model  Is  extremely  sensitive  to  the  PSD  slope 
of  the  terrain  profile  Input.  Of  significant  Interest  Is  the 
sensitivity  In  the  region  from  -2.0  to  -2.15,  since  many  natural 
terrains  fall  Into  this  region.  The  overlapping  of  the  lower  two 
lines  Is  probably  due  to  statistical  scatter  In  the  profile  generation 


process 
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II.  The  second  objective  of  this  research  was  to  study  the  sensitivity 
of  the  simulation  to  changes  in  profile  spacing.  One  random  profile 
was  created  using  NOIPSD  with  a  profile  spacing  of  0.5  inch.  One  may 
th  nk  of  this  profile  as  constant,  as  if  it  existed  In  Nature  somewhere. 
As  the  profile  spacing  (measurement  Interval)  Is  changed,  the  profile 
itself  will  not  be  altered.  What  will  change  Is  the  distance  between 
points  which  are  used  as  Input  to  the  vehicle  simulation.  The  first 
few  feet  are  shown  in  Figure  26(a).  Since  straight  lines  are  used  to 
Interpolate  between  point*  In  the  vehicle  simulation,  they  are  drawn 
between  them  in  the  figure.  If  a  one  Inch  interval  Is  desired,  every 
ether  point  is  "read"  by  the  computer  and  the  result  is  a  profile  as  In 
Figure  26  (b).  The  measurement  interval  Is  Increased  by  .5  Inch 
increments  In  Figures  26  (b-j).  Note  that  as  the  measurement  Interval 
i.  -reases,  the  profile  appears  to  become  smoother  as  the  valleys  are 
bridged  and  the  peaks  .  jttened  by  the  linear  interpolation  scheme. 

The  1/2- Inch- Interval  points  ore  repeated  in  Figure  26  (g)  for 
comparison.  Note  how  the  first  valley  has  been  gradually  disappearing 
and  the  peak  is  completely  flat. 

For  this  test,  an  RMS  elevation  of  1.0  Inch  was  chosen  In  order 
to  depict  a  more  natural  terrain.  A  PSD  of  slope  -1.850  was  specified 
to  keep  away  from  the  2.000-2.150  Instability  region.  The  desired 
profile  was  obtained  from  NOIPSD,  stipulating  NG  ■  7,  N  -  1800, 

RMS  -  1.0  and  y  ■  .046  In  the  input  phase.  The  final  PSD  equation  was 
PSD  -  .001134  (freq)"1,850 


(see  Figure  27). 


1  !■» 
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Each  vehicle  simulation  was  then  exercised  over  the  profile  at 
different  measurement  intervals.  Velocity  was  maintained  at  tan 
miles  per  hour.  The  simulation  was  exercised  until  a  steady-state 
absorbed  power  level  was  reached  and  the  average  absorbed  power  was 
used  to  compare  results.  Table  4  and  Figure  28  show  the  results. 

Note  that  the  average  absorbed  power  for  the  M35  maintains  an 
essentially  constant  level  throughout  the  tested  range  of  measurement 
Interval.  In  fact,  at  sixteen  inches,  the  tire  model  consists  of  only 
one  vertical  spring  —  a  "point  follower."  Thus  It  appears  that  for 
the  M35  truck  and  probably  for  other  heavy,  long-wheelbase  trucks, 
a  "point  follower"  tire  model  is  sufficient. 

For  the  MI51  Jeep,  however,  the  plot  shows  a  steady  decrease  in 
average  absorbed  power,  after  about  2.5  or  3.0  inches.  For  the  N 15 1 
and  probably  for  all  short-wheeibase,  small-tired  vehicles,  It  Is 
obvious  that  a  smaller  measurement  Interval  is  required  to  provide 
numerically  stable  absorbed  power  values.  The  data  Indicate  that  an 
Interval  of  as  low  as  2.5  inches  is  necessary.  For  profiles  with 
higher  RMS  elevations,  the  interval  may  be  even  smaller. 
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Table  4 

Average  Absorbed 
_ Power _ 


(watts) 


-HJiJL 

H35 

7.206 

1.328 

7.000 

1.207 

7.221 

1.296 

6.740 

1.276 

7.819 

1.377 

6.065 

1.234 

8.734 

1.318 

7.129 

1.404 

5.094 

1.254 

6.236 

1.186 

4.329 

1.082 

4.256 

1.016 

- 

1.321 
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CONCLUSIONS 

1.  The  vehicle  simulation  is  extremely  sensitive  to  PSD  slope, 
especially  In  the  neighborhood  of  -2.0. 

2.  The  vehicle  simulation  Is  sensitive  to  Input  profile  spacing 
for  light,  short  vehicles.  A  very  small  measurement  interval  of  from 
2.5  to  3  Inches  Is  necessary  for  numerically  stable  average  absorbed 
power  values. 

3.  The  vehicle  simulation  Is  relatively  insensitive  to  variations 
m  measurement  Interval  for  heavy,  large  vehicles.  Thus  a  much 
longer  measurement  Interval  is  possible  and  a  vastly  simpler  tire 
model  can  be  used. 
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RECOMMENDATIONS 

1.  That  tha  sensitivity  of  tha  modal  be  further  tested  using 
terrains  differing  In  PSD  slope  by  much  smaller  Increments  between 
-1.85  (or  even  -1.70)  and  -2.30. 

2.  That  the  M35  simulation  be  further  tested  for  sensitivity 
over  terrain  profiles  with  greater  AMS  elevations  and  different 
slopes. 
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APPENDIX  A(Documentatlon  of  NO  I  PSD) 

!.  DESCRIPTION:  This  program  Is  •  combination  of  Murphy's  two  programs 
N0I5E1  and  PSD.  It  first  generates,  using  Gaussian  noise,  a  random 
tarraln  profile  with  zero  mean  and  specified  RMS  elevation. 

II.  INPUTS  AND  OPERATING  INSTRUCTIONS: 

A.  Prior  to  execution: 

1.  The  equation  FLAG  -  6*N/(TAl>*62)  Is  written  assuming 

20  equivalent  degrees  of  freedom.  It  comes  from  Blackman 

4 

and  Tukey,  sections  B.23  and  B.24. 

2.  The  PSD  estimates  are  smoothed  by  hamming.  If  othar 
smoothing  Is  desired,  the  appropriate  parts  of  the  pro¬ 
gram  must  be  changed. 

B.  During  execution: 

1.  NG:  May  be  any  number,  but  to  Insure  differences  among 
profiles,  a  prime  number  should  be  used. 

2.  SHORT  TYPE  OUT:  Answer  YES  or  NO. 

a.  YES  causes  a  reduction  in  program  options,  setting 

TAU  ■  4.0,  RMS  -  4.0,  N»  1200,  and  no  output  files 
are  created. 

b.  NO  causes  the  program  to  run  normally,  with  all 

options  open. 

3.  In  put  of  TAU,  N,  and  RMS:  If  no  new  values  are  entered, 

the  program  will  set  them  to  4.0,  1200,  and  4.0 
respectively. 
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4.  ALPHA:  This  It  th«  partmtir  uttd  to  adjust  tho  slop* 
of  tho  PSD  curve. 

5.  FILES?  Antwor  YES  or  NO. 

o.  YES  eoutos  tho  program  to  ask  for  ft  la  nomas,  which 
ara  antarad  if  required.  If  no  antry  Is  mada  for 
althar,  that  ft  la  will  not  ba  opanad. 
b.  NO  causas  no  fllas  to  ba  craatad. 

6.  FLOW  •  :  Entar  tha  das  I  rad  low  cut-off  fraquancy.  If 
zero  Is  antarad,  tha  program  will  sat  FLOW  equal  to  .0002, 
which  normally  causes  only  tha  zero  fraquancy  estimate  to 
ba  deleted. 

III.  OUTPUTS: 

A.  PSD  file: 

1.  PSD  aquation  In  tha  form:  PSD  •  C(freq)n  , 

where  C  Is  tha  Intercept  at  fraquancy  •  1.0  on  tha  log- 
log  curvo.  With  PSD  on  tha  vertical  axis,  n  Is  tha  slope 
of  tha  curve  —  a  straight  line. 

2.  RMS:  This  Is  tha  actual  RMS  as  opposed  to  tha  desired  RMS. 

3.  Power:  Tha  area  under  tha  PSD  curve. 

4.  Number  of  Points:  Number  of  points  In  profile. 

5.  Cut-off  frequencies,  low  and  high. 

6.  Points  that  were  used  In  tha  square  fit  by  sequence  num¬ 
bers. 

7.  List  of  frequencies  and  PSD  estimates,  five  points  to  a 
line.  Only  tha  points  usei  in  tne  square  fit  ara  listed. 


8*  List  of  th«  points  not  ussd  In  ths  square  fit. 

8.  Road  flla:  Us to  as  Input  to  VEH. 

1.  Tha  profile  segment  spacing  In  Inches. 

2.  A  line  of  Identifying  information. 

3.  The  profile  points  In  order,  ten  points  per  line. 

C.  Teletype: 

1.  RMS. 

2.  PSD  equation.  (If  short  type  out  Is  selected,  only  the 
slope  of  the  PSD  curve  Is  printed  out;  none  of  the 
following  are  listed) 

3.  RMS,  power,  and  number  of  points  fitted  by  the  above 
equation. 

4.  SO.  (Standard  deviation  of  the  points  from  zero) 

5.  XBARAS:  Mean  of  the  profile  after  shift. 

6.  File  names  for  output  data.  (If  ---"-.DAT  is  typed  out, 
no  file  was  created.) 


IV.  SAMPLE  EXECUTION: 


A.  Normal: 


.execute  noipsd 

LOADING 

LOADER  9K  CORE 

execution 

NG  =  7 

SHORT  TYPE  OUT?  NO 

TAU=4.0,N=12.00,RMS=4.0  UNLESS  SPECIFIED  NOW* 

ALPHA  =  .0052 

LLAG  =  30 

FILES?  YES 

PSD  FILE  NAME:  PSD4 

ROAD  FILE  NAME*  R0AD4 

AUTO-COVAR.  IS  COMPUTED 

RMS  =  4.000000 

PSD  IS  FINISHED 

FLOW  =  .0002 

-2.000 

PSD  =  0.2067366E-02  FREW 

RMS  =  4.000000 

POWER  =  3.673574 

EQN  FITS  PTS  2  THRU  30 
'SD  =  0.6519456 

XBARAS  =  0.0000000 

FILE  FOR  PSD  PLOT*  PSD4  .DAT 
FILE  FOR  ROAD  PLOT*  ROAD4.DaT 


EXECUTION  TIME*  11.73  SE(  . 

TOTrtL  ELAPSED  TIME*  2  MI  N.  3.37  SEC. 
NO  EXECUTION  ERRORS  DETECTED 


B.  Abbreviated 


•EXECUTE  NOIPSO 
LOAOING 

LOADER  9K  CORE 
EXECUTION 

NQ  =  7 

SHORT  TYPE  OUT?  YES 
ALPHA  =  .0052 

SLOPE=-2.000 


EXECUTION  TIME*  6.40  SEC. 
TOTAL  ELAPSED  TIME*  24.75  SEC 
NO  EXECUTION  ERRORS  DETECTED 


PROGRAM  LISTING; 


inainwMi 


,iIMENS1GM  X(l2fltJ),DX(l2fc0)  .FREOf 
n I  HENS  ION  RX(25in,PX<25(9),SPX<?$0) 


WR|TE(6*90P 


ACCEPT  99/, NG 


•ft  '*■.  1  ) 


ACCENT  996,  OPTYPE 


'/RITE  (6,902) 


LL AG8 IP  I Y ( PL AG+ • 5  ) 


FORMAT ( * ♦LI  AG  8  '  »  I  A , /  ) 


FORMAT ( '♦FILES?  »,S) 


tF(OPFlL,EO, 'NO' )  GO  TO  190 


FORMAT! ’♦PSD  FILE  NAME  I  »,S) 


IF(OPTYPE,EQ, ’YES’ )  GO  TO  190 


4081 

rw7t*  » 

500 

n  *•*< 

»•  SJGMAN  AND  ST,  DEV.  ARE  COMPUTED 

41  c* 

700 

SIGMAN*RMS*S0RTM..-EXP(-2t»AI  PMA«TAU)) 

520 

AA«EXP(-ALPHA*TAU) 

538 

SD*SIGMAN*SIPMAN 

540 

C  *•*< 

»•  PISPLACEMEMTS  ARE  OBTAINED 

558 

c  •  ••< 

»•  FROM  GAUSSIAN  RANDOM  NUMBERS 

560 

SUM*0, 

•5781 

nX(N)a0 . 

5B0 

V  N  *  N  “  1 

*590 

'.IJMlstf. 

600 

no  210  1*1. NN 

A15L_ 

CALL  GAURNU(V.NG) 

*27 

nX< 1 )*V 

630 

210 

SUMsEUM*DX( I ) 

640 

XBAR  =  StlM/NN 

650 

HO  270  1*1. NN 

660 

nxm*nx(l  >-xbar 

670 

220 

suMi*suMi»nx< i > 

6BP 

xbaras*sumi/nn 

690 

eact*i. 

700 

Y«0. 

710 

X  ( 1 )  *  0  . 

72? 

no  230  1*2, N 

730 

230 

X ( 1  )  *DX ( I > ♦X  < I -1 ) *A A 

740 

C  •••< 

►i  RMS  IS  COMPUTED  AND  ADJUSTED 

750 

C  *••< 

••  TO  DESIREO  LEVEL 

760 

240 

SUMX<*0 , 

770 

00  250  1*1. N 

780 

xm*FACT»X(I) 

790 

250 

SUMX?*SUMX2*X( I  )  *X  ( I  ) 

_Bidl2 _ 

ARMSsSORT  <SUMX?/'N) 

010 

EaCT*RmS/AHmS 

020 

Y*Y*1, 

830 

IE'Y.EO.1.)  GO  TO  240 

040  C  •••( 

»•  AUTO-COVAHI ance  computation 

850  260 

MP0 1 N  T  *  N 

860 

no  2*0  I*1.LLAG 

870 

$X*0 , 

6B0 

M*I-1 

890 

f.  Xs.NPOINT-M 

9B0 

.  FNX*rLQAT (NX) 

910  DO  27(8  J*1,NX 

921 _ K*M  +  J _ 


930  270  SX*SX*X( J)«XCK> 

_ Z8J?. _ EX ( I >  fSX/FHX 


950 

l  t  .  r  x  m  "  ...  - . 

RMSiSCRT (RX<1 )  ) 

960 

IF (DPTYPE ,nE . ’YES *  > 

WRITE (6,912)  RMS 

970 

91?  FOR"»T< ’♦auto-covar. 

'  RMS  =  ' »  G10 . 7/ ) 

IS  COMPUTEO ' / » IX » 

I 

I 

I 
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i 

* 


i 


* 

i 


t 

i 


• « 

!> 

•  «• 

ii 

ii 

ii 

ii 

ii 

fi 


194? _ me*\f*i 


195? 

196? 

924 

WR I TE ( 21  *  924 ) 

FORMAT (/// *  THE  FOLLOWING  POINTS  WERE' 

197? 

♦ 

'  not  FITTEOJ •/) 

1913? 

IF(NR.EO.l)  GO  TO  41? 

199? 

WRITE (21, 922)  ( < FREQ < K > » SPX ( K ) ) , K«1 ,NS*1> 

2??? 

r-0  TO  42? 

2?1? 

41? 

WRITE (21, 925) 

2?2? 

_  925 

FORMAT!’  NO  LOWER  FREQUENCIES  WERE  DELETED*/) 

2?3? 

42? 

IF (NF ,EQ,LLAG>  GO  TO  460 

2040 

WRITE (21, 922)  !(FREQ(K).SPX(K)).KsMr.LLAG) 

205? 

GO  TO  470 

206? 

463 

WRI TF (?1 , 926 ) 

2?7? 

926 

FORMAT!/’  NO  HIGHER  FREQUENCIES 

WERE  DELETED’//) 

??e? 

.4  7? 

CONTINUE 

2?9? 

C 

21  fc!? 

48? 

I F  !  F  N 1 .  E  o ,  •  *)  FW1*’ - - 

211? 

IF  (FM2.E0*  ’  *)  Fn-2»* . » 

212? 

IF(OPTYPE.NE.’YES’)  WRITE(6.927)  SD. 

213? 

i 

XBARAS.FM1.FN2 

214? 

927 

FORMAT!’  SO  *• #F12*7*/’  X8ARAS  *’. 

2150 

♦ 

F12 » 7 . / *  FILE  FOR  PSD  PLOTI  ’ .  A5 . ’ , 0 A T ’ , / 

216? 

♦ 

’  FILE  FOR  ROAO  PLOT  *  • . A5. * .DAT ’ / > 

217? 

!F!Fvl.EO, • . . .  Go  TO  49? 

218? 

WRI  TE  (21. 920  )  ARMS»SD.X0AR,XBARAS.‘,G.N 

219? 

928 

F0RMAT(67X. ’RHS  *’*Fi2.7,/ 

220? 

♦ 

53X. ’STANDARD  DEVIATION  *’#F12.7,/ 

'221? 

♦ 

54X.  ’  X-BAH  AFTER  GAUSS  *’,F12.7,/ 

'222? 

♦ 

?4X.’X-UAR  AFTER  SHIFT  «*,Fl2.7*/ 

223? 

♦ 

45X, ’START ING  NUMBER  FOR  RAN ( 2 )  *’*I4,/ 

224? 

♦ 

55X  ,  ’  NUMBER  OF  POINTS  «’,I6) 

225? 

929 

FORMAT(1?F12.0) 

226?__ 

49? 

TF(F* 2.E0. • . . .  GO  TO  9999 

227? 

WR I TE ( 22  *  9  30  >  T  AU.  FN2,  AM,  NG,  ALPHA,  T  All.  RMS 

228? 

93? 

F0RMAT!1X,F/1X.A5,’ , DAT- -SLOPE* ’ » F6 , 3, 

229? 

♦ 

•  MGs’.Ic*’  ALPHA*’ ,F6. 5, 

230? 

♦ 

♦  T AU* ’ * E 3 . 1 » ’  RMS**.F3,1) 

231? 

no  5?0  l*m* m, 1? 

232? 

50? 

WRITE (22*929)  (X (K) * Ks | -9, 1  ) 

233? 

9999 

CALL  EXIT 

234? 

996 

E0RmAT(A5) 

235? 

997 

FORMAT! I ) 

236? 

998 

FORMAT(E) 

237? 

999 

FORMAT!/) 

138* 

END 

239? 

2400 

R.produc.d  from 
betl  .v.il.b!#  copy. 
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APPENDIX  B( Documentation  for  VEH) 

I.  DESCRIPTION:  This  program,  in  its  original  form,  was  also  obtained 

from  N.  R.  Murphy,  WES.  It  Is  a  non-linear,  time-domain  wheeled 
vehicle  simulation  on  a  digital  computer,  a  POP  10. 

II.  INPUTS  AND  OPERATING  INSTRUCTIONS: 

A.  Prior  to  execution:  The  vehicle  constants  must  be  entered 

in  Subroutine  DATA  (see  Appendix  C)  and  the  tire  load- 
deflection  curves  must  be  on  hand. 

B.  During  execution: 

1.  Name  of  input  file  (Road  file  created  with  NOIPSD). 

2.  Desired  Delta-L;  This  should  be  a  multiple  of  the 
Input  profile  spacing. 

3.  Name  of  vehicle  to  be  simulated. 

4.  Options.  Answer  YES  or  NO. 

a.  Absorbed  Power. 

b.  Detailed  output  file? 

(1)  YES  causes  program  to  ask  for  file  name,  which 

should  not  be  more  than  five  characters  long. 

j 

(2)  NO  causes  program  to  skip  the  ab^ve  question 

and  no  output  file  will  be  created. 

c.  Peak  accelerations. 

d.  Driver  motions  (only  asked  if  vehicle  is  a  track) 

e.  RMS  of  all  accelerations.  (This  must  be  answered 

YES  If  PWRPIT  is  to  be  used  later  to  plot  the 

i 

! 

i 

i 
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absorbed  power. 


f.  Tire  deflections.  The  program  will  type  out  the  loads 


on  each  axle,  front  to  rear.  Enter  the  deflection 


caused  by  each  load  on  a  single  tire. 


5.  Vehicle  velocity  In  miles  per  hour. 


6.  Teletype  printout  time  Interval  In  seconds. 


7.  Time  and  absorbed  power  only?  Anser  YES  or  NO. 


a.  YES  causes  the  program  to  type  out  only  time,  absorbed 


power,  and  average  absorbed  power  for  each  Interval. 


b.  NO  will  cause  the  normal  printout  format  to  be  used. 


8.  Stop?  Answer  YES  or  NO.  Program  will  stop  on  entry  of 


III.  OUTPUTS:  These  depend  on  the  options  selected.  Assuming  all 


are  desired,  the  outputs  will  be  as  follows: 


A.  To  detai led  file: 


I.  List  of  THRESH  and  GAM'IA. 


2.  Vehicle  velocity  In  MPH  and  inches  per  second. 


3.  Delta-L. 


4.  Delta-t, 


5.  Number  of  steps  in  RKG  integration. 


6.  RKG  step  size. 


7.  The  input  profile  I de  ltl f (cation. 


8.  Time,  profile  input  point,  average  power,  absorbed  power, 


and  vehicle  motions  at  each  step. 
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6.  To  teletype  Beslcally  the  seme  information  as  above  Is  print¬ 
ed  out  on  the  teletype,  except  that  the  time,  Input  point,  etc. 
are  printed  out  only  according  to  the  teletype  prlnt-out  inter¬ 
val  specified  earlier. 

IV.  SAMPLE  EXECUTION: 
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IV.  SAMPLE  EXECUTION: 
A.  Norma  1 : 


.EXECUTE  VEH 
FORTRAN!  VEH 
LOADING 

LOADER  9K  COKE 
EXECUTION 

FILE  NAME  OF  INPUT  PROFILE!  ROADX 
DESIRED  DELTAL:  1. 

NAME  OF  VEHICLE?  M151 

DO  YOU  WANT  THE  FOLLOWING  OPTIONS? 

ABSORBED  POWER?  YES 
A  DETAILED  OUTPUT  FILE?  YES- 
FILE  NAME!  FI  LEX 
PEAK  ACCELERATIONS?  NO 
RMS  OF  ALL  ACCELS?  YES 

FROM  FORCE-DEFL  CURVES  FOR 
M— 151  JEEP  TIRES,  ENTER  DEFLECTIONS* 
LOAD  :  581.071  .82 

LOAD  =  623. 249  .865 


VEHICLE  VELOCITY  IN  MPH*  10.0 
TTY  PRINTOUT  TIME  INTERVAL!  .2 

TIME  A  POWER  TYPEOUT  ONLY?  NO 


VELOCITY: 10.00  MPH  (176.30  IPS) 

DELTA-L:  1  .030  DELT.l-1 :0.0057 
NS T EPS:  5  H:  .001  136 

VEHICLE  IS!  M-I51  JEEP 
INPUT  PROFILE  IFi 

ROADX.DAT— FILE  FOR  VARIABLE  MEhSUREMENT  INTERVALS 


1 


! 

* 


v 


ft- 106 
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A.  Normal:  (continued) 


DISPL 

VELOC 

mCCEL 

RNSAC 

TIME= 

0.000  INPUT: 

:  0.000 

ABSORBED 

POWER: 

<;-G 

>1.17601 

0.00000 

0.00000 

0.00000 

PITCH 

>0.00086 

0.00000 

0,00000 

0.00000 

AXLE1 

-0.82060 

0.00000 

0.00000 

0.00000 

AXLE2 

-0.86500 

0.00000 

0.00000 

0.00000 

STOP? 

NO 

TIME= 

0.205  INPUTS 

-0.209 

ABSORBED 

POWER:  l 

C-G 

-1.22561  - 

3.25482 

-0.21239 

0.15562 

PITCH 

•0.00026  • 

0.13723 

-2.62326 

2.13225 

AXLEI 

-1.0692R  - 

8.41071 

1.23916 

0.45973 

AXLE2 

-0.86037 

1.79340 

0.04520 

0.09609 

STOP? 

NO 

TIMES 

0.403  INPUTS 

-0.549 

ABSORBED 

POWER:  1 

c-u 

-1.09423  - 

1.51579 

0.12303 

0.23038 

PITCH 

0.00089  > 

0.05403 

2.84386 

3.34899 

AXLEI 

-0.72493  - 

0.32713 

0.13655 

0.65754 

AXLE2 

-0.80245 

0.23447 

-0.15274 

0.12382 

STOP  7 

YES 

EXECUTION  TIMEi 

29. 

27  SEC. 

0.140  0.041 


0.515  0.231 


TOTAL  ELAPSED  TIME*  3  MIN.  57.30  SEC. 
NO  EXECUTION  ERRORS  DETECTED 


EXIT 


R-1  Gk6 
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B.  Abbreviated: 


.EXECUTE  VEH 
LOADING 

LOADER  9K  CORE 
EXECUTION 

FILE  NAME  OF  INPUT  PROFILE*  ROADX 
DESIRED  DELTALi  1.5 
NhME  OF  VEHICLE?  M151 

DO  YOU  WANT  THE  FOLLOWING  OPTIONS? 

ADSORBED  POWER?  YES 
A  DETAILED  OUTPUT  FILE?  NO 
PEAK  ACCELERATIONS?  NO 
RMS  OK  ALL  ACCELS?  NO 

FROM  FORCE-DEFL  CURVES  FOR 
M— 151  JEEP  TIRES,  ENTER  DEFLECTIONS* 
LOAD  =  581  .071  .82 

LOAD  =  <523.249  .86r> 


VEHICLE  VELOCITY  IN  MPH*  10.0 
TTY  PRINTOUT  T I I1E  INTERVAL*  .1 


TIME 

&  POWER 

TYPEOUT 

ONLY?  YES 

TIME 

ABSPWR 

AVEPWR 

0.00 

0.00 

0.00 

0.10 

0.05 

0.01 

0.20 

0.19 

0.04 

0.31 

0.73 

0.18 

0.40 

0.78 

0.30 

0.50 

0.61 

0.38 

0.61 

0.88 

0.43 

0.71 

0.75 

0.50 

0.80 

0.89 

0.52 

0.90 

1.G4 

0.64 

STOP? 

YES 

EXECUTION  TIME* 

39.78  SEC. 

TOTAL 

ELAPSED 

TIIC* 

1  MIN.  57.87  SEC. 

NO  EXECUTION  ERRORS  DETECTED 
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V.  PROGRAM  LISTING: 


STORAGE  ALLOCATIONS  I 

NOTE  I  THE  FOLLOWING  COMMON  STATEMENTS 
MUST  APPEAR  IN  EACH  SUBROUTINE 

- fcOMHON  F6RCHTAr.FORCT(7).rORClCtAl.rCRCyC61 _ 

COMMON  SPDEF (6) »0SPSF (6) # THRESH (200) « SIGMA (9)  ! 

- COMMON  VARMA)»Y(260)iPWRVAR<9)«OAMP(6) _ 

COMMON  ACC!$S(9>  # ACCGSI9) #ACCMAX(9) » ACCMIN'9) 
COMMON  SUMRMS(Q).RMS(9),LEN(lfl),HASS(61 - 

COMMON  h,t.qeltat>oeltal>velips.velhph,nsteps 
COMMON  Y INI.  QPVMAX.DRVMIN.ABSPMR, GAMMA (200) 
COMMON  0 1 SORV , VELORV , ACCORV , RHSDffV 

- COMMON  1FPWR»  IFFILE,  1FPACC,  IFDRWa.1.F  RMB - Sf^l 

COMMON  NY, I  OF  *  NAXLES, NSEGS* IFhORZ.FNAME 


COMMON  VEHQIO(2^PROFIL(26«)»PASTPC260)  *  INDEX  ! 


v 


COMMON  SUM  JT  (4,3) ,SSL0PF(5,3),SINT(5»3) 
"common  01  1MIT(P.3).DSL0PE (3,3), DINT (3i3J 
niHENSlQN  0RIVER(4)>I0PT(6),NYTEHP(3I 
DIMENSION  FI0(12) *XTNAME(4) 

DIMENSION  FK(ia),P(lB).0(l8)tPV(9),pwRrK 

DIMENSION  PP(9) ,QQ(9) 

RE*L  LEWiijAffg,  IWRT1A _ 

EQUIVALENCE  USETUPiNY) 

EQUIVALENCE  CPWIVEWC11 .OISQRV) _ 

EQUIVALENCE  ( I  OPT (l)i JFPWR ) 

nATA  XTNAME/»M1S1». t,»  «/ _ 

NO»»N» 

IYES«»Y» _ 

IFSTOP«*N* 

mmiL  initialization _ 

DO  100  I >1 « 9 

PWRVARI I)«g. _ 

00(1 )*0, 

ACC15S(I)«0. _ 

sumrms  < i >  *0 , 

JUiilBj _ 

ACCMAX ( I ) *0 , 

ACCMINC  I  »«B, _ . _ 

continue 

ALCDRV«0. _ 

DO  110  I *1 i 18 

VAP(I)«0, _ 

T»0, 

SOVRMS»0. 

AB$PWR«0, 

DRVMAX*0, 

DRVMIN*0 • 

Y 1 N«fc , 


I 


R-1 6h6 


85 


500 

510 

H*.001 

0ELTAL»4, 

'520 

151EL 

JST0**1 

VSTOP«0 

P  ''Lj 

1540 

1550 

C  MM. 

TPR1MT«0, 

OPEN  INPUT  FILE  AND  READ  PROFILE 

1560 

'570 

C  M«M 

SPACING  ANb  OTHER  IDENTIFICATION 
WRITH6.900) 

•560 

590 

ForhaT('  >ILE  name  of  lNPU'T“PR6riLL'i  ',!) 
REA0I5.996)  FN2 

600 

610 

IF(FM2.EQi'  •>  GO  TO  120 

CALL  IFILEI22.FN2) 

620 

630 

61fl 

c  ••••• 

120 

CALL  FlLIfTTFID* JS .SPACING .1)  ’ 

VEHICLE  CONSTANTS  READ-IN 

WRITH6.902) 

650 

660 

902 

FORMAT! ’♦NAME  OF  VEHICLE?  f#$) 

REA0(5»  996 )  TNAME 

670 

660 

00  130  I>1.4 

IF!TNAME-XTNAME!I))1 30 *140*130 

690 

130 

continue 

WR  T  TE  ( ft  < 903 )  XTNAME 

710 

720 

903 

* 

FORMAT!’  THE  AVAILABLE  VEHICLES  ARE | * 

4 ! 2Xi A5 ) ) 

730 

740 

140 

GO  TO  120 

GO  TO  !150< 16 0*170*180)  ! 

750 

760 

150 

ASSIGN  400  TO  I SUB 

GO  TO  190 

770 

760 

160 

ASSIGN  410  TO  ISUB 

GO  TO  190 

790 

600 

170 

ASSIGN  420  TO  ISUB 

GO  TO  190 

010 

620 

160 

190 

ASSIGN  430  TO  ISUB 

CONTINUE 

830 

040 

C 

SELECTION  OF  OPTIONS 

HR  I TF (6  »  904 ) 

650 

660 

904 

« 

FORMAT!'  DO  YOU  WANT  THE  FOLLOWING'* 

'  OPTIONS?'//'  ABSORBED  POWER?  '.!) 

870 

660 

REA0<5,995)  IFPWR 

WRITE16.905) 

690 

960 

905 

FORMAT ! ' *A  OETAILED  OUTPUT  FILE?  ’iS) 

READ(5i 995)  IFFILE 

910 

9?fl 

IF !  1FFILE.EQ, 'N' )  GO  TO  210 

WR ! TE ! ft »  906  > 

930 

910 

906 

FORMAT! '♦  FILE  NAMEl  »,S) 

ACCEPT  996.  FN1 

950 

960 

210 

CALL  0FILE!21,FN1) 

WRITE(6»907) 

970 

sag. 

907 

FORMAT! '♦PEAK  ACCELERATIONS?  ',S) 

RE ADI 5. 995)  IFPACC _ 

IF  ( 1  *LT  .  3 )  CO  TO  215 


♦90 


[V.T  . 


FORMAT! '♦DRIVER  MOTIONS?  '.$> 


GO  TO  216 


1050  216  WR I TF ( 6#  909 ) 

1060  909  FORMAT  (  '  ♦RMS  OF  ALL  ACCELSt  f«S) 

1070  R£AD(5, 995 )  IFRM* 

0 _ 


C  «*••#  INITIALIZE  PROFILE  ARRAYS 


PASTP  < I >  *0 


PROFILd  )*0 


NYT£Mp(2)»aEN(i)^LEN<2)-LEN(3))/OELTAL 
NYTrMP(3)a(LFN(n*LEN(?)  +  LFN(4 


VEHICLE  RUN  VARIABLE  INPUT 
WRI TE  < 6 #  911 > 

F 0 R M A T ( / ♦  VEHICLE  VELOC 


REA0I5.998)  VELMPH 
HR  I  TE ( 6  #  912 } 


FORMAT ( * ♦TTY  PRINTOUT  TIME  INTERVAL!  ».$> 
REAO(5i 998 )  TIP 


WR J  TE (6# 


REA0(5#995>  OPTYPE 

TIME  step  AND  rkg  time  SET-UP 


VEL I pS= VELMPH* 17 • 6 
AT«DELTAL/VELIPS 


NSTEPS«DELTAT/H 

TEMp.NSTEPS 


H«OELTAT/TEMP 


nO  2<0  1*1* IDF 
ACCGSn  )«ACCISS<  I  )/38 


ACCGS(2>*ACC1SSC?> 


ABSPWR=0 


CALL  AVERACE(ABSPWR#AVEPWR) 
FdFORV.EQ.  »N'  >  ' 


0ISDRV«VAR(1)^DRVLEN»VAR(2J 


"M5  CALCULATION 
F(lFRMS,EQ. «M» ) 


CO  TO  280 


R-1646 


147? 

1A&£ 

1490 

im l 

151? 
152  L 
1530 
154? 
155? 
1560 
1570 
1560 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670" 
1660 
1690 
1700 
1710 

JLZ2JEL 

1730 

1740 

1750 

1760 

1770 

1780. 

1790 

1600 

lfllE 

'1620 


PO  270  l«l.IDF 

SUHRMS(1)«SUHRHS 

0HSm«0, 

ir(T,Ne.e>)  rhs< 

IFUFDRV.EO,*N») 

SDVRMS«SDVRM5»ACI 

RMSORV«0* 


GO  TO  300 


PEAK  ACCELERATION  CALCULATION 
IFUFPACC.EQ,  »N«)  GO  TO  320 

no  310  1*1 f IDF 

ACCMAX(I)«AHAX1<ACCHAX(I )«ACCCS(I) ) 
ACCMlN(I)>AMTNl(ACCM!N(I)'ACCGSn)) 
IFUFQRV.EQ.  «N* )  GO  TO  340 

0RVMAX*AMAX1(DRVMA*» ACCDRV) 
DRVKIN»AHIN1(0RVM1N»ACC0RV) 

PROFILE  INPUt 

CALL  FIIIN(FI0,JST0P#SPAC I NG . 2) 


C  PROGRAM  OUTPU1 


IFdFFlLE.EQ.  «YM  CALL  FILWRTIFIO.NPL. 

"Fn^avepurj -  — - 

f F < T • L T iTPRINT )  GO  TO  360 

CALL  PRINTIFID.  trsW.OPTYPE.AVEPWR) 

TPRINT«TPKlNT*TtP 

tFMFSTOP,F6,  *  Y 1 )  GO  tO  510 

MAIN  PROGRAM 

CALL  SHIFT 

SHIFT  ADVANCES  THE  Y  PROFILE  ARRAY 
IF (J$toP.LQ,2)  GO  Td  500 
PROFtL<l)«YlM 

■InCiEX—KSTEFS - 

LUF»?»IDF _ 

DO  380  J»l# IDF 
K* J* I  OF 

PY(J)»VAR<K) - 

HO  390  I*1»NY 

V( I)»PASTP( I )♦( < INDEX*NSTEPS*114 _ 

(PROFIL(I)-PASTPC I )  M/NSTEPS 

no  440  I « 1 , 4  _ 

GO  TO  !SUB» (400,410*420,430) 

CALL  WHEELS (FKiNVTEMP) - 

GO  TO  440 


♦  11* 


» _ _ . _ _ _ . _ _ _ 


'  * 


u 


R-I&6 


100  _11? 


180  130 


SUBROUTINE  FILIN <F  ID, JS.SPAC INC, N> 
THIS  SUBROUTINE  READS  A  NEW  INPUT 


DIMENSION  FYIN<10> 

.DATA  EllMltiSuZ _ 

IFIN.LT.2)  GO  TO  130 

-IY»  I  Y*MM _ 

IFUY.GT.10)  GO  TO  120 

Y.INf.Fmmi _ 

return 

-RE  A  0,122x900)  (FYINlDxli 
FORMAT  < IFF ) 

..If  IEQPCI..GQ  TO  999, g _ 

I  y«  i  Y-ie 


URITE(6,901) 

FORMAT ( 1 ♦  DESIRED  S 
RE AD<  5 , 998 )  OELTAL 
,-BEAQI  22x30  2.1  ..SPA  Cl  J 
( F Y I N I J ) , J*1 # 10 ) 

FOPM*T  lllli  F/JLXiA2  A! 
mmi-deltal/spacing 

JJY _ 

RETURN 

J5.*2 _ 

RETURN 


SUBROUTINE  RUNGE(P.0»X,FK,H,N) 
c  THIS  SUBROUTINE  IS  THE  RKG  ALGORITHM, 

- DIMENSION  P(1?,Q(1?.X(1),FK(1),A(4),B(4),C(< 

DATA  A/, 5», 292893219, 1 , 70710678, ,166666667/ 

_ HAJA-l/jxil,flM,,2.«^ _ 

DATA  C/. 5, ,292693219. 1.70710678,, 5/ 

_ TA«A(N) _ _ _ 

TB*8<N> 

_ TCiC(N) _ ~ 

DO  100  I *1 » M 

_ E.Ln«TAMFKQ.LrTB>QUn _ 

x ( i »»x ( I >*p  c  i  > 

JJEUa _ Qm»QIl)*3.»P(l)>TC«rK(II _ 

RETURN 


R- |6U6 


jMMMLt 


SUBROUTINE  OATA(NI) 

THIS  SUBROUTINE  CONTAINS  THE  VEHICLE 

Parameters;  it  also  CALLS  GAMSUB - 

DIMENSION  I  SETUP ( 5>  ).i()Sl  GMA<.?i  *JL 

nt  HEMS  I  ON  0VEHCL(P.4)  .nTHRSH(9.4).DGAMHA(J 

DIMENSION  OHASS(6,4) ,0VAR(9,4) ,DLEN(10»4) 

n I  HENS  1  DM  OFHASSt4).DlNRTA14)JDORtfLNlAJ _ 

01  MENS  I  ON  USL1M(4,3,2),DSSL0(5,3,2) 

DIMENSION  OS  I  NT (5. 3,2) _ 

DIMENSION  00L1M(?#3,2> ,DOSLO(3 ,3,2) 

DIMENSION  DDINT(3.3,2) - 

DIMENSION  OR ( 2 ) 

INTECEP  DSETUP_(5._4I - 

REAL  LENiMASS » I  NR  T I A 


DATA  DVEHCL/10HM-151  JEEP,10HM-35  TRUCK, 

1  ~  1  Pi u m[- 6 0  TANK  i.lflHfl  _ _ 

OSETUP'S  ARE  NY,!DF,NAXLES,NSEGS,IFH0R2 

- DATA  Q£ETUP/g,4.7,0,fl.g.5.3.fl.fl. _ 

t  50, 9, 6, 5, 1,36, 8,  r>,  5,1/ 

_ DATA  DTHRSH/9*0..9«0.. _ _ 

1  3 , 5 , 1  •  ,  0  , ,  1  •  , 3 , 5 » 4 *0 . , 


DATA  DGAMMA/9*0.,9«0., 

-1 - 3aaS.,4  7l5J-..SflgB..47lS..3HS^..4>a.. _ 

1  1500, ,2000. >350 0. ,2000. ,1500. >4 *fl,/ 

_ OATA  DSIGMA/9*0.  .9*0.  . _ _ 

1  3145, ,1670  ,0, ,-1670. , -3145 ,, 4*0 , , 

J - I5flfl.i7fl0«>fl.«-7y»a.  .4»W./ _ _  _ 

DATA  0LEN/44.3>4C!.7,8«a,,H3.>39,,24f  ,24,  ,6«|0., 

1  77.. 44., 11. ,-22., -55. 4*0., _ 

1  52.,24.,0,#-28»>"65.,5»0,/ 

DATA  DhASS/.27i.27i _ 

1  1.191,2.08,2.05,3*0. >6*0, ,6*0./ 

_ OAT  A  -QE.MASS/2.5flil8.g.iZ.ig»/ _ 

DATA  0 1 NRT A/3282. ,90876. ,0, ,0./ 

_ DATA  DORVLN/0,i0.,25.,25./ _ 

DATA  DVAR/9*0, ,9*0. , 

1  -5. 79. -.D089. -.966. -.97. _ 

1  **  •  942,  -  .  913,  -  ,  884 ,  —  ,  856, 0  ,  , 

1  »3.75, -.0087. -.76.-.78.".76,-.73,-t6fi«2*0./ 

data  0?L  I M/12*999, ,-4.4,-3.65,3.65,4.4, 

1  -5. 7. -5.1 1 5.1 1 5, 7, -5. 7, -5.1. 5, 1,5.7/ _ 

DATA  DSSLO/15#1500. , 

1  11771 .43, 3333.33.1145.2,3333. 33, 11771.43, 

1  46000., 9333, 33, 2509, 8,9333.33, 46000,, 


OATa  OSiNI  715*0,  , 

1  66714, 4t 7986. 65. 0.0,-7986.68,-66714,; 

1  2438*0. # 34800, ,0, 0.-34800, .-243800 i* 

1  243800., 34800., 0.0, -34800., -243800,/ 

OATA  DDL  I  M/4 *999, #2*0 1 » 

OATA  DOSLO/9*42»  , 

1  70. .1402. .40. .0. .1 963.. 0. .8. .1883.10 

OATA  OOINT /9*0 . » 

1  »6M.  «0.  »620.  »»958.  »8.  #951*  >«98l«  >8.«j 

OATA  DR/15.,19,/ 

DO  170  1*1,2 _ _ _ 

VCMQfOC  1 )  ■OVEHCLd  »N> 

DO  110  1*1,5 _ _ _ 

!  SETUP  (  I)  ■OSETUPU  ,N)  - 

DO  120  I *1 , IDE 

_ LEN<n«DLENtliN) _ 

VAR  C !  )«OVAR{J ,N) 

_ 00  130  1 *1 ,NAXLES _ 

MASS(I)>OHASS(!,N) 

_ FHASSiDFHASS(N) _ 

INRT1A*D!NRTACN> 

_ ORVLEN«PPRVLM<N) _ 

IF(N.LE,2)  GO  TO  150 

_ 00  140  1*1, MSECS _ 

THRESH  <  1  )*OTHRSHU  »N) 

_ S1GMA(I)*PSIGHA<1.N) _ 

GAMHA(I)«OGAMHA(!,N) 

_ RETURN _ 

00  200  1*1, NAXIES 

•  00  160  J»1.4 _ 

SLIMITtJ# I)«OSLIH(J,  I.N) 

_ DO...  170  J«li? _ 

SSLOPE<J,I)*OSSLO(J»I,N> 

_ SINT(J,I)«DSINUJi.I»Nj _ 

00  160  J*l,2 

_ nLlHlT(J,n»00LlH(J,ltN) _ 

DO  190  J*t,3 

_ DSLQPE ( J « I ?,DDSLQ£ J/ li  N? _ 

OINT(J,I)*ODINT<J,I,N> 

_ CONIINUS _ 

R«ORCN> 

_ CALL,  GAHSUB(N.R) _ 

DO  210  l*l,NAXLES 

_ VAR.(  I*2)*~YY ( I ) _ 

WEIGHTC I > ■WEIGHT( I )-MASS ( 1 )*386, 

_ IF (NAXLES «LT.3)  GO  TO  210 _ 

IFC.LT.2)  GO  TO  210 

_ wElQHTU?«2.*HElPHT<n _ 

SPOEFU  )«-HEIGHT(I  )/SSL0PE(3, 1 ) 


"  RWWpn  m 


T1«VAR(3)»SP0EF(1) 


T2«<T2*LEN(4)*T3«LEN(3) )/(LEN(3)*LEN(4>) 


VAR<2)«(T2-T1)/WB 


8 

IP 

FORMAT (5F10  *5) 

RETURN 

R- 16*6 


93 


die 

020 


SUBROUTINE  GAMSUB<N,H, 

PI  MANSION  VEHWT (2) #  RF ( 2  > i WDF ( 2 )  * WQR ( ?#  2) 


nJl" 

|040 

1050 

iwmp/gw.i _ 

DIMENSION  DELTA<0/90> 

REAL  LEN.MASS 

1060 

070 

DO  100  1*0*30*1 

D(M*I*OELTAL 

080 

>090 

100 

110 

ir<<0<I)/R>.GT, 0,788)  GO  TO  110 

VSEGS* ! *2"1 

•  100 

NY*(LEN<1)*LEN(2)*LEN(4) )/OELTAL*NSEGS 

iiir 

KK* I -1 

•120 

KZEROsI 

0sLFN(l)«»LEN(2) 

140 

150 

UE!GHT(1)«(FMASS«(LEN(2)/0)4MASS(1))*3R6, 
WEIGHT (2 ) * ( ( (FMASS* (LEN(  1 1 /Q J )/ (NAXLES-1) ) 

160 

<70 

6 

♦MASS  <  2 ) ) *386 , 

TF (MAXiESiLT.3)  GO  TO  120 

180 

1  90 

WEIGHT ( 2 ) *WE I GHT (2 )/2 , 

WEIGHT  (3  )*(<(FMASS*<LEN(l)/Q))/<NAXLES-in 

200 

210 

6 

120 

+MASSI3) )»386,/2, 

WRITE<6i900)  VEHOID 

220 

230 

900 

♦ 

FORMAT!  '  FROM  FORCE-DEFL  CURVES  FO«  '/IX. 

2A 5.'  TIRES,  ENTER  DEFLECTIONSl  •/) 

240 

950 

00  135  I*1iNAXLES 

WR I TE ( 6 , 901 )  WEIGHT ( 1 ) 

260 

27* 

901 

135 

FORMAT (  ’ ♦  LOAD  *  * , FB . 3 , JX , * ) 

RE AO ( 5 , 902 )  YY < I ) 

'280 

i290 

902 

FORMAT (F ) 

SUM*is0. 

•380 

*310 

SUM6*0 . 

SUM7*0. 

'320 

330 

DO  130  I *0#KK  » 1 

nPRIMEU  )*SQRT(R**2.-D(I>**2.  > 

340 

350 

COSTHETA*DPRIMEU)/R 

THRH  > -R-OPRIME < I ) 

360 

•370 

0ELTA<M*YY(1)-THRU> 

IF (DELTA ( I ) ,LT  ,0,0)  DELTA(I)*0.0 

360 

SUM5sSUM5»OELTA(I)*COSTHETA 

390 

J*  I  ♦•<«♦  1 

400 

410... 

0ELTA(J)«YY(2)-THR(I ) 

IF(DELTA(J) .LT,0,0)  D£LTA{ J) *0,0 

420 

SUM6*SUM6+OELTA(J)*COSTHETA 

430 

JJ*J*KK*1 

440 

450 

0ELTA(JJ>*YY(3)-THR<  I  ) 

IF(OELTA(JJ).LT.C,0)  DELTA (JJ ) *0,0 

460 

470 

137 

SUM7*S')M7  +  0ElTA(JJ)*C0STHETA 

CONTINUE 

460 

CSTSIS 

490 

SUM5»SUM5*2 

,-D 

ELTA(0) 

500 

SUH6*SUM6*2 

.-OELTA(KK*i) 

csssts 


csssis 


SPKF«WEIGHT(1)/SUM5 

SPKR1«WEIGHT(2)/SUM6 


SPKR2*0,0 

F(N.EO.l)  GO  TO  136 


SPKR1«2*SPKR1 
SPKR?«2* WEIGHT  (3)/Sl)M7 


00  140  I*KK,0,-1 


THRESH(K?ERO-I )*THR<  I ) 
THRESH(KZ£RQ*I ) «  T»  R  jjj 


DO  i50  I *1 » NSEGS 

THRESH (  I  ♦NSEGS >* THRESH < I ) 


THRESH  <  J ) *THRESH (I ) 


HO  160  I*KK,P,-1 


JJ»KHERO*I 


GAMMA< J)»SPKF»COSTHETA 


Jl* J*NSEGS 


GAMMA ( Jl)«SPKRl*COSTHETA 


J2*J1*MSECS 


GAMMA ( J2)»SPKR2«C0STHETA 


890  .9 

OHfll  9491 

32 

910  9 

922 _ 

33 

FORMATI&PC  D(  ' » I2» ' >  *  »F8.2/)/) 


FORMAT  ( // '  I  THRESHU)’ 


00  240  I «1 * 2*NSEGS 


32  FORMAT (1X»I2#1X#2F15»3) 


33  F0RMAT(/6X, »REAR| • ) 


WRITF(?1#934) 


00  260  I«2*NSEGS*1.3«NSEGS.l 


RETURN 


4 


SUBROUTINE  PEAKAC(NPL) 
THIS  SUBROUTINE  WRITES 


•MAXIMUM  MINUMUM*) 


FORMAT ( *  C-C  *2F9,4,/»  PITCH  S2F9.4) 


IBKFVIH! 


FORMAT < '  AXLfi Ii#lX»2F9,4) 


♦  ACCMINUOF) 


IFMFORV.EO.'Y*)  WRITE<0.905)  DRVMAX.ORVMlN 


IFUFFltE.EQ.’N')  RETURN 


URITE<21i9B2)  (ACCMAXU)iACCMIN(!)iI«1i2) 


IFUFHORZ.EO.il  WRITE(21»904)  ACCMAX(IDF)* 


lF{|FORV»EQt  * V f )  HRITE(21#905)  ORVMAX.ORVMfN 


SUpROUTINE  AVERAGE *X, Y ) 

THIS  SUBROUTINE  COMPUTES  THE 


OATA  N/0/ 


II  1JI  Jlh 


Y«SUM/N 


E 


97 


460 

*oa 

190 

IF<IFDRV,EO.*N»>  60  TO  160 

uftfrr<A.909)  fmnvnm.fai.ifi 

900 

009 

FORMAT ( 1  O*!VC*'|4F10,9) 

UL. 

160 

HRITE(6.910i 

920 

910 

FORMAT ( *  6T0R?  *»S> 

«0 

ACCEPT  996.  t FSTfiP 

940 

996 

FORMAT(Al) 

«« 

_KK«0 

960 

RETURN 

171 

170 

_IFlfriRST.0Q.0l  URITE(6.9111 

i960 

911 

FORMAT<//'  TIME  ABSPHR  AVEPWR*) 

399 

IFIRSTal 

400 

WR!TE(6i912)  T, ABSpHR. AVEPHR 

A10 

012 

FORMAT (1X.F9.2.2F7.2) 

620 

KK"KK*1 

630 

IFfKK.CE. 101  00  TO  160 

440 

END 

010 

SUBROUTINE  VARFIX(I.VOUT) 

020 

C  •  *••* 

THIS  SUBROUTINE  IS  CALLED  BY  PRINT 

aaa 

C  •••it 

TO  SELECT  THE  VARIABLES  TO  RE  PRINTED. 

040 

DIMENSION  VOUT (4 ) 

•00 

VOUT ( 1 ) sVAR  f  f  1 

060 

N» I ♦ I  OF 

070 

VOUT ( 2 ) *VAR ( N) 

060 

VOUT ( 3 ) ■ACCCS ( I ) 

000 

VOUT  t 4 ) «RMS l I 1 

100 

RETURN 

Ilf  , 

_EM2 - 

*-)6k6 


I 


98 


I 

I 

0 

0 

0 

0 

0 


0 

0 

& 


*10 

020 

gat 

C 

SUBROUTINE  POWER(FK) 

THIS  SUBROUTINE  CALCULATES  ABSORBED  POWER* 
PIHENSIQN  FK<9> 

040 

050 

*  ^  f  . 

U2««67'743*ACCDRV-1'042*PWRVAR(6) 

U1«-U2-3.246*PWRVAR(6) 

060 

070 

U0«-Ul4l,318*PWRVAR<4> 

FKfl)«H»< .00873*PWRVAR(2>*PWRVARf3>> 

080 

690 

FK<2)4HM**4 , 9948 4* ACCQR V) 
FKf3>«H*t-100.*U&-59.*PWRVARf3>) 

100 

410 

♦ 

FK(4)«H»<-13.*U1*71,6*PWRVAR(5> 

-53.49*PWRVAR(4)> 

120 

<30 

FK(5>*H*(-100.*U1-47,780U0> 

FKf6>«H»<-10.»U2-78.59«PWRVAR(7> 

140 

150 

♦ 

-55,28*PWRVAR(6>) 

FKt7J«H»<-l0.*U2-6.259*Uii 

160 

470 

ifab*  -  , 

■  4 

FK<8>»H«<-677,43*ACCDRV*388.8*PWRVAR(9> 

-46.67*PWRVAR<B>> 

160 

190 

FK(9>«H«(-67.743»ACCDRV-2,742*U2) 

RETURN 

200 

,r 

END 

4JL0L 

_SliflROUTlNE  SHIET 

200 

21SL 

c  •«••• 
_ c  MMt 

THIS  SUBROUTINE  ADVANCES  THE  PROFILE 

UNDER  THE  VEHICLE  *  ^ 

220 

23JL 

DO  100  I«NY,2,-1 ! "  - - 

■PASTPU-1  )«PRoni  <i-i) 

240 

•2S1- 

100 

PR0FIL(I)«PR0FILU-1)  - - - 

-RETURN 

260 

END  - - - 

0 

& 

0 

0 

0 

fl 


R- 1 61*6 


010  SUBROUTINE  FILwRT(nD,NPL»FN2»AP) 

020  c  •••••  THIS  SUBROUTINE  HANDLES  THE  OUTPUT 


030 

C  ***•« 

TO  AN  EXTERNAL  FILE 

040 

OIMENSION  HEADK6) » HE AD2 ( 2 > 

050 

DIMENSION  VOUT (10)  #F1D(12) 

060 

DATA  IFIRST/0/ 

070 

DATA  HE ADI/ ' AXLEl  ’  #  ’ AXLE2’ # ’ AXLE3 • # ' AXLE4  ’ , 

060 

♦ 

•AXLES’# ’AXLE6’/ 

090 

DATA  HEA02/’H,C-G’.’V.ORV’/ 

100 

IF ( IF IRST  #  NE • 0  >  GO  TO  i30 

110 

WR I TE  <  21 # 902 )  VELMPH# VELIPS.DELTAL.DELTAT. 

120 

* 

NSTEPS.H 

130 

902 

FORMAT!/#’  VELOCITY  *  ’.F6.2# 

140 

♦ 

•  MPH  <  ’  F6 . 2 » ’  IPS)  ’  ,7X, ’OELTA-L*’ 

<50 

♦ 

*F5.3* »  INCHES, * .OX* ’OELTA-T*'. F10.B. 

160 

♦ 

•  SECONOS’#//’  RUNGE-KUTTA-GILL  INTEGRATION* 

170 

* 

NUMBER  QF_  STEPS*’. I4#10X. ’STEP  SIZE  (H)*’# 

160 

♦ 

E12.6//) 

<90 

WRITEI21.903)  VEHOID.EID 

200 

903 

FORMAT!//’  VEHICLE  IS:  ',2A5#// 

210 

♦ 

*  PROFILE  IS*  ' * 12A5 ) 

220 

IFIRST  *1 

230 

KK*2 

240 

IF ( IFORViNE* 'N*  >  GO  TO  120 

250 

HEA02 (2 ) *0  > 

260 

KK«KK*1 

270 

120 

IF ( IFHOR? • NE • 0  * )  GO  TO  160 

260 

HEA02(1)*HEAD2(2) 

290 

KK«KK-1 

300 

GO  TO  160 

310 

130 

IF (NPL.LT ,50)  GO  TO  170 

320 

140 

IF (NPL*CE # 54 )  GO  TO  150 

33P 

VR! TE (21#  904 ) 

340 

904 

FORMAT  < 1H  ) 

350 

NPLsNPL*l 

360 

GO  TO  140 

370 

150 

NPL*P 

360 

160 

WRITE (21# 905)  (HEAD1(I),I*1,NAXLES)» 

390 

♦ 

(HEA02 ( I ) # I *1 , KK ) 

400 

905 

FORMAT ( 1H1 » '  TIME  Y ( I ) ' • 14X# 

410 

♦ 

♦V.C-G  PITCH’ #4X.8(A5.4X)) 

420 

WRITE(21#904) 

430 

NPL«NPL*2 

440 

170 

DO  180  I »1 # IDF 

450 

180 

VOUT ( I )*VAR  ( I  ) 

460 

J«!OF 

470 

IF(IFORV.EO.’N’)  GO  TO  190 

460 

J*J*1 

490 

VOUT(_U)»i)ISORV 

R- 16U6 


101 


010 

020 

C  #***/ 

SUBROUTINE  WHEELS ( FK  #  NYTFM0 ) 

THIS  SUBROUTINE  IS  THE  GENERALI 2FD 

01301 

yum  rn  vehirle  morel 

040 

050 

DIMENSION  TEMPC3) ,FK(i8),NVTEMP(3) 

REAL  len.mass# inrtia 

060 

070 

e 

OATA  IFIRST/0/ 

RESULTING  AXLE  FORCES 

080 

090 

99 

DO  110  I*t.NAXLES 

F  QPCW ( 1 ) s0 , 

100 

110 

00  1O0  Jsl.NSECS 

Jjs J+NYTEMP ( I ) 

120 

*30 

JJJsJ*(I-1)«NSEGS 

T£Mp0=Y< JJ)-VAR(2+I ) -THRESH (JJJ) 

140 

150 

100 

IF  (TEMP0.LT, P. )  TEHP0S0.P 

FOPCWd  )*FORCW(  I  ) ♦GAMMA  ( JJJ )#TFMP0 

160 

110 

CONTINUE 

170 

_ C  M«»« 

180 

190 

UU*LFN ( 3  )  ♦LEN  <  4  ) 

IF(UU.LE.0.)  GO  to  120 

200 

210 

U*(VAR(4)-VAR(5) )/UU 

RETA*ATAM(U) 

220 

230 

08ETAs(VAR(4*I0E  >-VAR(5*IDE))/{UU»(l.*U*U) ) 
TEMP?sSIM<9ETA) 

240 

?50 

120 

TEMP4=C0S<BLTA) 

TEMpl*SI N ( VAR ( 2 )  ) 

260 

970 

. . 

TEMP3*C0i <VAR(2> ) 

SUSPENSION  SPRING  DEFLECTION  fSPDEF) 

280 

290 

SP0EF(1)*VAR(3)-VARC1>-LEN(1)*TEMP1 

SP0EF(2)*VAR(4)-VAR<1)*LEN<2)*TEMP1 

300 

IS10 

8 

-LEN ( 3 ) *TEMP2 

IF (NAXLES iLT.3)  GO  TO  130 

<320 

330 

4 

SPOEF {3)*VAR(5)-VAR(1)  ♦LENI2 )*TEMP1 
♦LEN(4)*TEHP2 

340 

350 

c  •• 

130 

SUSPENSION  SPRING  RATE  OF  DEFLECTION  (DSPDF 
0SPDF<l)*VAR(3MDF>-VAR(i*IDF> 

360 

370 

8 

-LEN<1)*VAH2*IUF)»TEMP3 

nSP0F(2)*VAR(4+IDF)-VAR(i*I0F) 

380 

390 

8 

♦  LtN<2)»VAR(2*I0n*TEMP3-LEN(3)*0BFTA*TEMP4 
inNAXLES.LT. 3)  GO  TO  140 

400 

410 

OSPDF(3)eVAR(5*IOF)-VAR<l*IOF> 

♦LEN(2)*VAR(?*IDF)*TEMP3*LEN(4)«RBETA*TEMP4 

420 

420 

C  *• 
140 

•  •• 

SUSPENSION  SPRING  FORCES 

00  160  I «1 » NAXLES 

440 

450 

K  =  5 

DO  150  M»l»4 

460 
4?g . 

!F<SPDEHI>,LT.SLIMIT(M,I)>  GO  TO  155 
CONTINUE 

460 


GO  TO  160 


:  . 
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C  •• 
c  *• 


FORCK<n«SSLOPE(K#I)«»SPDEFei>*SINT<K#n 

SUSPENSION  DAMPING  FORCES _ 

00  181  i-l.NAXLES 

K»S _ _ _ 

00  170  H«1.P 

|FCDSPOFC|),LTtOHHITCH,m  GO  TO  17S 

continue 

GO  TO  180 _ 

K*M  ~  " 

,nAHP(I)«OSlOPE(K»I)»DSPOF<I)*PINTCK>n 
DIFFERENTIAL  EQUATIONS 

FK 1 1 )  AND  FK(1*IDF)“-»CG  MOTION _ 

FX<2>  ANO  FK«2*IDF>-  — PITCH  MOTION 
FK<3>  ANO  FK(3*I0F>  — AXLEl  MOTION 

FK<N)  ANO  Fk!n*IDF>— AXLEN  MOTION 


FKU>«H«VARU*IOF> 

stemp»p. _ 

00  2*0  I ■! « NAXLES 
_ lEMili  "F.Q5.CK  t 1  i  tOAHP  ( I ) _ 

STEMP»STEMP*TEMPU> 

_ FK.M*2»1,0Q«H«(F0RCH(I)-TEHI 

t  “MASS (I )*386 t ) /MASS ( I ) 

_ FKllMPFlSHgtSIEaP-rfilASStMi 

FK(2*I0F)«H*(LEN(1)*TEMP(1) 
.1 _ sLE.Ht2l.#TEMP(2}  I/INRTIA  ,  . 

RETURN 


THIS  SUBROUTINE  IS  FOR  A  TRACKED 


SUBROUTINE  M60(FK) 

IMIS.  SUBROUTINE  li-Lfl&JU 

VEHICLE*  NOT  USEO  HERE 


It 
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APPENDIX  C: 

The  vehicle  parameters  necessary  for  wheeled  vehicle  simulation 
ara  listed  In  the  same  order  as  in  subroutine  data.  The  variable 
names  are  written  in  all  capital  lettars.  See  Figures  Cl  through  C3 
for  clarification.  Those  parameters  In  DATA,  not  listed  herein, 
are  not  necessary  for  wheeled  vehicle  simulation.  This  appendix  Is 
intended  as  a  work  sheet  for  use  In  obtaining  data  for  simulation  of 
a  new  vehicle. 

A.  DSETUP 

1.  NY  -  number  of  profile  points 

under  vehicle  (computed)  .  0.0 

2.  IDF  -  degrees  of  freedom . . 

■  4  for  2-axle  vehicles 

■  5  for  3-axle  vehicles 

3.  NAXLES  -  number  of  axles  ...  .  . 


4.  NSEGS  -  number  of  segments  In  tire 

model  (computed) . .  0.0 

5.  IFH0RZ . 0.0 

B.  DLEN 

1.  LEN(l)  -  horizontal  distance  from  front 

axle  to  center  of  gravity . ___  Inches 

2.  LEN(2)  -  i2:  horizontal  distance  from  CG 

to  rear  axle  (or  center  of  rear 

bogie  assembly) . .  _ __  Inches 

3.  LEN(3)  -  £g:  horizontal  distance  from 

center  of  rear  bogie  assembly 

to  center  of  2nd  axle  (■  0  for 

2  axles)  .  inches 
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4.  LEN(4)  -  it  horizontal  distance  from  center 
of  rear  bogle  to  center  of  3rd  axle 


(■  0  for  2-axles) . .  _____  Inches 

5.  LEN(5)  thru  LEN(9)  .  0.0  Inches 

C.  DMASS  -  Unsprung  masses 

1.  MASS(l)  -  H.:  Hass  of  front  axle  In-  lb-sec* 

eluding  wheels . . inches 


2.  MASS(2)  -  Mgj  Hass  of  second  axle 


lb-sec* 

Inches 


3.  HASS(3)  -  Hj:  Hass  of  third  axle 

4.  HASS(4)  -  HASS(6)  . 

0.  DFHASS-Ho:  Sprung  mass  . 


Ib-sec* 

inches 

Ib-sec* 

inches 

Ib-sec* 
I  nc'hes 


E.  DINRTA-I  j  Pitch  moment  of  Inertia  of 
spring  mass  about  CG  .  .  . 


Ib-sec* 

Inches 


F.  DESCRIPTION  OF  SUSPENSION  SPRING  FORCE 
FUNCTION  (see  Figure  C2) 

I.  DSLIH  -  limits  of  regions 


a.  Front  axle: 

(1)  S  L I  HI  T(  1,1)  -  x,, . Inches 

(2)  S  L I  HI  T(  2, 1)  -  Xg, . Inches 

(3)  SLIHIT(  3, 1 )  -  x3J . Inches 

(4)  SLI HIT( 4, l)  -  Inches 

b.  Second  axle: 

(1)  SLIHIT(I,2)  -  x,2 . Inches 

(2)  SLIHIT(2,2)  -  Xgg . Inches 


fitrea  (lbs) 


Segmenttd  llntar 
approximation 


Equations  of  lines  are: 

F,j  -  mjjdj  ♦  C,j 

I  ■  region  number  (»l,2, 3,*»,5) 
J  «  axle  number  (1,2,3) 


Figure  C2.  Segmented-Linear  Spring 
Force  Approximation 


H-l  646 


(3)  SLIHIT(3,2)  -  Xjg . lncl»«* 

(4)  SUHIT(3,2)  -  x^ . Inchaa 

c.  Third  axis  (If  2-axlas,  that*  ara  xeroea) 

(1)  SLIHIT(  1,3)  -  *13 . Inchaa 

(2)  SLIMIT(2,3)  -  Xg3 . Inchaa 

(3)  SLIHIT(3,3)  -  x33 . Inchaa 

(4)  SLIHIT(4,3)  -  Inchaa 


2.  OSSLO  -  Slope  of  Itnea  which  approximate  the  force 
va.  deflection  curve  of  the  suapenaton  aprtnga. 

a.  Front  axle: 

(1)  SSLOPE(  1,1)  -  m,, . . 

(2)  SSL0PE(2,  l)  -  i^| . . 

(3)  SSL0PE(3,I)  -  . . 

(4)  SSLOPE(M)  -  "V! . . 

(5)  SSL0PE(5,1)  -  . . 

b.  Second  axle: 

(1)  SSLOPE(  1,2)  -  m,2 . . 

(2)  SSL0PE(2,2)  -  n^g . . 

(3)  SSL0PE(3,2)  -  . . 

(4)  SS10PE(4,2)  -  . . 

(5)  SSL0PE(5,2)  -  . . 

c.  Third  axle: 

(1)  SSLOPE(  1,3)  -  «"13 . . 

(2)  SSL0PE(2,3)  -  1^3 . . 

(3)  SSLOPE(3,3)  -  ^ . . 


nKgin 

' 
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(V)  SSL0PE(4,3)  -  "V3  •  •  . _ 

(5)  SSL0PE(5,3)  -  «^3 . _ 

3.  DSIHT  -  Intercepts  of  fore*  vs  deflection  lines. 


Front  exle: 


(1)  SINT( lf I)  -  cn 

(2)  SINT(2,|)  -  Cg, 

(3)  SINT(3,1)  -  c31 

(4)  sint(4,i)  - 

(5)  SINT(5,1)  -  c5, 


b.  Second  exle: 


(1)  SINT(I,2)  -  C|g 

(2)  SIMT(2,2)  -  Cgg 

(3)  SIMT(3,2)  -  c^ 

(4)  SINT<4,2)  - 

(5)  SINT'  2'  -  c„ 


Third  t 


(1)  SINT(1,3)  -  cu 

(2)  SINT(2,3)  -  «23 

(3)  SINT(3,3)  -  c33 

(4)  SIMT(4,3)  -  c^3 

(5)  SINT(5,3)  -  c53 


6.  DESCRIPTION  OF  SUSPENSION  DAMPING  FORCE  FUNCTION 
(see  Figure  C3) 


4'  v  :  •  , 

.  »  ^  V 


I.  00LIM  -  Limits  of  regions 


c.  Front  axle: 

(1)  DUMIT(  1,1)  -  vn . 

(2)  DLIMIT(2, 1)  -  v2|  . 

b.  Second  axle: 

(1)  DUMIT(I,2)  -  vJ2 . 

(2)  DUMIT(2,2)  -  . 

c.  Thl rd  axle: 

(1)  0LIHIT(  1,3)  -  v13 . 

(2)  0LIHIT(2,3)  -  v23  . 

2.  DDSLO  -  Slope  of  lines  which  approximate  the 
force-velocity  curve, 
e.  Front  axle: 

(1)  DSL0PE(1,I)  -  n,,  . 

(2)  DSL0PE(2, 1)  -  ng,  . 

(3)  0SL0PE(3,1)  -  n3, . , 

b.  Second  axle: 

(1)  DSL0PE( 1,2)  -  n|2  . 

(2)  DSL0PE(2,2)  -  . 

(3)  0SL0PE(3,2)  -  n^ . 

c.  Third  axle: 

(1)  DSL0PE(  1,3)  -  nJ3 . 

(2)  0SL0PE(2,3)  -  n23 . 

(3)  0SL0PE(3,3)  -  . 


3.  OOINT  -  Intercepts  of  linos  on  deeping  force 


vs*  deflection  velocity  curvet, 
e*  Front  exle: 

(1)  OINK i,i)  -  dn  .  .  .  . 

(2)  0INT(2,  l)  -  . 

(3)  0IMT(3*  I)  -  d3,  .  .  .  . 
b*  Second  exle: 

(1)  DI*T(I,2)  -  d|2  .  .  .  . 

(2)  DINT(2,2)  -  *  *  .  * 

(3)  0IMT(3,2)  -  d^  .  .  .  . 

c.  Third  exle: 

(1)  DIHT(I,3)  -  dJ3  .  .  .  . 

(2)  0IMT(2,3)  -  d23  .  .  .  . 

(3)  0INT(3,3)  -  d33  .  .  .  . 


H.  DR  -  The  undeflected  tire  red  I  us  In  Inches 


Inches 


APPENDIX  0  (Documentation  of  PWRPIT) 

I*  DESCRIPTION:  This  program  plots  tha  absorbed  power  end  the 
average  absorbed  power  against  time.  It  stores  the  plot  In 
e  file  on  disc  which  can  be  printed  on  a  line-printer  with 
130  characters  per  line. 

II.  INPUTS  AND  OPERATING  INSTRUCTIONS: 

A.  Prior  to  execution:  none 

B.  During  execution: 

1.  I  FILE:  The  name  of  the  file  which  contains  the  power 
and  absorbed  power.  This  file  must  be  In  the  format 

O. 

of  the  detailed  output  file  from  VEH,  with  options 
absorbed  power  and  RMS  accelerations  specified.  This 
file  name  must  be  5  or  tass  characters  long. 

2.  OH  IE:  Desired  name  of  plot  fIIej  not  over  $  characters 

III.  OUTPUTS:  One  file  containing  a  plot  of  absorbed  power  and 
average  power  vs.  time. 


SAMPLE  EXECUTION: 


.EX  PVRPLT 
FORTRAN*  PVRPLT 
LOADING 

LOADER  AX  COKE 
EXECUTION 

IFlLEl  FILEX 
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OFILEi  PWRX 

END  OF  FILE  ON  DSK1 

C//F,/8X,Ft/7X,F) 

LAST  FORTRAN  1-0  AT  USER  LOC  000267 


EXECUTION  Tl.'lEl  B.20  SEC. 
TOTAL  ELAPSED  TIME*  45.20  SEC. 
NO  EXECUTION  ERRORS  DETECTED 


IF  ( I XX  •  CT « 11(1  >  |XX*11I 


WR!TC(22)930)  TU  ) ,  <X<  J) ,  , 


FORMAT ( * *' »F5»IfH0Al#2F6tf ) 


XC!XX)«BUNK 


WRITC(22»949)  (!,!>», 22) 


